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Dredged  material  is  often  disposed  Into  diked  containment  areas  (DMCAs).  Topically, 
the  dredged  material  is  pumped  into  a  containment  area  as  a  slurrv  that  is  ID-  to 
2D-percent  solids  by  weight.  Water  and  solids  separate  in  the  containment  area  because  of 
gravity  sedimentation,  and  the  clarified  water  is  discharged  over  a  wc- 1  r  and  through  a 
culvert  to  a  receiving  water  body.  TTie  suspended  solids  concentration  of  DMCA  effluent  Is 
Inversely  related  to  the  residence  time  of  water  in  the  DMCA.  The  longer  the  water  stavs 
in  the  DMCA  before  It  is  discharged,  the  lower  the  effluent  turbidity  and  suspended  solids 
concentrat ion  will  be.  Residence  time  is  therefore  an  Important  consideration  in  design  of 
DMCAs  . 

Designers  attempt  to  sire  DMCAs  to  provide  adequate  sediment  storage  and  eftluent 
quality  at  least  cost.  The  <(irrentlv  recommended  |iroredure  for  siclng  I'Mi  As  uses  the 
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19.  ABSTRACT  (Continued). 


estimated  average  residence  time.  Prior  to  this  study,  estimates  of  mean  residence  time 
were  based  on  only  a  few  observations  of  operating  DMCAs.  The  mean  residence  time  was 
simply  estimated  to  be  0.44  times  the  ponded  volume  divided  by  the  inflow  rate. 


This  study  presents  observed  residence  time  distribution  data  from  some  12  DMCAs. 
These  data  were  collected  by  Injecting  Rhodamine  WT  dye  at  the  inflow  point  and  monitor¬ 
ing  dye  concentration  in  the  outflow.  The  resultant  time-concentration  curve  is  identi¬ 
cal  to  the  residence  time  distribution  during  the  dye  test.  In  addition,  similar  data 
for  chlorine  contact  chambers,  physical  models,  and  waste  stabilization  ponds  found  in 
the  literature  are  also  examined. 


A  composite  data 
the  ratio  of  mean  res 
time  is  equal  to  the 
width  ratio.  Baffles 
ratio  of  a  given  DMCA 
overall  size.  The  re 
incorporated  into  the 
shape,  and  spur  dike 


set  from  all  different  sources  showed  a  strong  relationship  between 
Idence  rime  to  theoretical  residence  time  (theoretical  residence 
pond  volume  divided  by  the  average  flow  rate)  and  the  pond  length- 
or  spur  dikes  can  be  used  to  increase  the  effective  length-to-wldth 
,  thereby  increasing  the  mean  residence  time  without  increasing  the 
latlonship  between  mean  residence  time  and  length-width  ratio  was 
DMCA  design  procedure,  and  a  process  for  selecting  DMCA  size, 
layout  to  achieve  minimum  cost  resulted. 


The  utility  of  state-of-the-art  mathematical  hydrodynamic  models  for  simulating  flow 
through  DMCAs  is  examined  in  Appendix  A.  It  was  found  that  the  mathematical  models  are 
capable  of  reproducing  the  observed  residence  time  distributions  fairly  well  if  model 
coefficients  are  arbitrarily  adjusted  to  achieve  best  fit.  However,  better  field  data 
will  be  required  to  develop  reliable  a  priori  modeling  capability.  Although  additional 
research  based  on  higher  quality  field  data  is  needed  to  provide  more  definitive  informa¬ 
tion  on  DMCA  residence  time  distribution,  at  this  point  it  Is  apparent  that  in  many  cases 
DMCA  economics  could  be  Improved  by  the  use  of  baffles  or  spur  dikes. 
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PREFACE 


This  study  was  conducted  as  part  of  the  Dredging  Operations  Technical 
Support  (DOTS)  Program  at  the  US  Army  Engineer  Waterways  Experiment  Sta¬ 
tion  (WES),  Vicksburg,  Miss.  The  DOTS  Program  is  sponsored  by  the  Office, 

Chief  of  Engineers,  US  Army,  through  the  Dredging  Division  of  the  Water 
Resources  Support  Center,  Fort  Belvoir,  Va.  The  DOTS  is  managed  by  the  WES 
Environmental  Laboratory  (EL)  through  the  Environmental  Effects  of  Dredging 
Programs  (EEDP) .  Dr.  Robert  Engler  was  Program  Manager  for  the  EEDP,  and 
Mr.  Thomas  R.  Patin  was  the  DOTS  Program  Coordinator,  Mr.  David  B.  Mathis  was 
the  WRSC-D  Technical  Monitor  during  preparation  and  publication  of  this  report. 

This  report  was  prepared  between  April  1985  and  April  1986.  Data  analy¬ 
sis  and  preparation  of  the  main  body  of  the  report  were  done  by  Mr.  F.  Douglas 


Shields,  Jr.,  Dr.  Edward  L.  Thackston,  and  Dr.  Paul  R.  Schroeder  of  the  Water 


Resources  Engineering  Group  (WREG) ,  Environmental  Engineering  Division  (EED) , 
EL.  Technical  support  was  provided  by  Mses.  Cindy  Cox,  Kathy  Smart,  and  Anita 
Zitta,  all  of  the  WREG.  The  DYECON  computer  program  mentioned  herein  was 
written  by  Mr.  Steve  Pranger  of  the  WREG,  and  Mr.  Pranger  assisted  with  use  of 
the  program  for  this  study.  Many  helpful  suggestions  were  made  by  Dr.  Thomas 


Walskl,  also  of  the  WREG.  Dr.  Walski  and  Dr.  Michael  R.  Palermo,  Chief,  WREG, 


acted  as  technical  reviewers  for  the  draft  report.  WREG  personnel  worked 
under  the  direct  supervision  of  Dr.  Palermo  and  Mr.  Clifford  L.  Truitt,  Acting 
Chief,  WREG,  and  under  the  general  supervision  of  Dr.  Raymond  L.  Montgomery, 
Chief,  EED,  and  Dr.  John  Harrison,  Chief,  EL. 

Field  data  used  in  this  study  were  collected  by  WREG  personnel  and  Brian 
J.  Gallagher  and  Company,  contractor.  References  cited  herein  may  be 


consulted  for  the  names  of  those  involved.  Wind  data  were  obtained  from  the 


National  Weather  Service  of  the  National  Oceanic  and  Atmospheric  Administra¬ 


tion  except  as  otherwise  noted. 

The  study  described  in  Appendix  A  was  performed  by  Messrs.  Donald  P. 


Bach,  William  D.  Martin,  and  Samuel  B.  Heltzel  of  the  Estuarine  Simulation 


Branch,  Estuaries  Division  of  the  WES  Hydraulics  Laboratory  (HL)  .  Mr.  Robert 


A.  Boland,  Jr.,  was  Chief  of  the  Estuarine  Simulation  Branch.  General  super¬ 


vision  was  provided  by  Mr.  William  H.  McAnally,  Jr.,  Chief  of  the  Estuaries 
Division,  and  by  Messrs,  H.  B.  Simmons  and  Frank  A.  Herrmann,  Jr.,  who  were 
Chiefs  of  the  HL  during  the  study  period.  Appendix  B  was  prepared  by 
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Dr.  Schroeder  of  the  WREG.  This  report  was  edited  by  Ms.  Lee  T.  Byrne  of  the 
WES  Information  Products  Division. 

COL  Allen  F.  Grum,  USA,  was  the  previous  Director  of  WES.  COL  Dwayne  G. 
Lee,  CE,  is  the  present  Commander  and  Director.  Dr.  Robert  W.  Whalin  is 
Technical  Director. 

This  report  should  be  cited  as  follows: 

Shields,  F.  D,  et  al.  1987.  "Design  and  Management  of  Dredged 
Material  Containment  Areas  to  Improve  Hydraulic  Performance," 

Technical  Report  D-87-2,  US  Army  Engineer  Waterways  Experiment 
Station,  Vicksburg,  Miss. 
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CONVERSION  FACTORS,  NON-SI  (METRIC) 

UNITS  OF  MEASUREMENT 

Non-SI  units  of  measurement  used  in  this  report  can  be  converted  to  SI 
(metric)  units  as  follows: 


Multipl 


To  Obtain 


acres 

4,046,873 

square  metres 

acre-feet 

1,233.489 

cubic  metres 

cubic  feet 

0.02831685 

cubic  metres 

cubic  yards 

0.7645549 

cubic  metres 

degrees  (angle) 

0.01745329 

radians 

feet 

0.3048 

metres 

gallons 

3.785412 

cubic  decimetres 

Inches 

2.54 

centimetres 

miles  (US  statute) 

1 .609347 

kilometres 

square  feet 

0.09290304 

square  metres 

square  yards 

0.8361274 

square  metres 

yards 

0,9144 

metres 
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DESIGN  AND  MANAGEMENT  OF  DREDGED  MATERIAL 
CONTAINMENT  AREAS  TO  IMPROVE 


HYDRAULIC  PERFORMANCE 

PART  I :  INTRODUCTION 


Background 


1.  The  US  Army  Corps  of  Engineers  TCE)  annually  supervises  the  dredging 

3* 

of  about  480  million  yd  of  sediment  from  the  Nation's  rivers  and  harbors 
(Francingues  et  al.  1985).  Due  to  environmental,  economic,  or  Institutional 
reasons,  a  large  quantity  (though  fractionally  small)  of  this  dredged  material 
must  be  placed  on  land  in  diked  dredged  material  containment  areas  (DMCAs) . 

In  DMCAs,  solids  in  the  dredged  material  are  separated  from  water  by  sedimen¬ 
tation,  and  the  resulting  clarified  water  may  be  discharged  back  to  the  water¬ 
way.  DMCAs  typically  have  irregular  shapes,  range  from  10  to  100  acres  in 
area,  and  have  ponded  depths  of  2  to  10  ft.  Inflows  to  these  areas  are  from 
dredge  pipes  and  typically  range  from  5  to  30  cfs.  Inflows  tend  to  be 
unsteady  and  Intermittent.  Outlets  from  DMCAs  are  usually  weirs  with  stop-log 
crests.  The  hydraulic  performance  of  a  DMCA  is  reflected  by  the  frequency 
distribution  of  residence  time  for  water  passing  through.  The  estimation  of 
residence  time  is  an  important  factor  in  the  design  of  DMCAs. 

2.  Preliminary  guidance  for  determining  residence  time  was  developed  by 
Montgomery  (1978)  during  the  Dredged  Material  Research  Program  (DMRP) . 
Refinement  and  validation  of  the  procedures  were  conducted  under  the  Disposal 
Operations  Technical  Support  (DOTS)  program.  These  efforts  included  field 
studies  of  residence  time  in  a  number  of  DMCAs.  This  report  describes  the 
results  of  the  field  studies  and  the  development  of  refined  techniques  for 
estimating  residence  time  of  DMCAs.  The  effects  of  various  design  decisions 
and  management  techniques  on  residence  time  distributions  are  also  explored. 


Use  of  Residence  Time  in 
Design  of  Containment  Areas 


3.  Montgomery  (1978)  developed  a  design  procedure  for  DMCAs  for  fine¬ 
grained  dredged  material.  This  work  showed  that  design  must  be  based  on  two 
general  phenomena:  (a)  the  settling  characteristics  of  the  dredged  material 
and  (b)  the  hydraulic  characteristics  (residence  time  distribution)  of  flow 
through  the  DMCA. 

4.  The  settling  characteristics  of  dredged  material  are  a  function  of 
many  variables  and  must  be  experimentally  determined  using  laboratory  proce¬ 
dures  developed  by  Montgomery  (1978)  and  extended  by  Palermo  (1984).  Whether 
the  slurry  undergoes  flocculent  settling  (the  manner  clays  normally  settle  in 
fresh  water)  or  zone  settling  with  a  distinct  Interface  (the  manner  sediments 
normally  settle  in  salt  water),  the  clarif ication  (removal  of  suspended  solids 
by  settling)  is  a  strong  function  of  time.  A  typical  set  of  curves  describing 
flocculent  settling  of  a  slurry  at  different  initial  suspended  solids  concen¬ 
trations  as  a  function  of  time  is  shown  in  Figure  1.  Obviously,  the  longer 
the  residence  time,  the  more  solids  are  removed  by  settling. 


SS„  INITIAL  SUSPENDID  SOLIDS 
CONCENTRATION 

SS  -  SS  •  SS 


Figure  1  .  Removal  of  suspended  soIiiL-  with  (  irri- 


5.  The  ideal  type  of  residence  time  distribution  in  a  DMCA  would  be 
plug  flow.  Plug  flow  may  he  described  as  a  situation  in  which  each  parcel  of 
water  entering  the  basin  proceeds  through  the  basin  in  exact  order  of  entry, 
occupies  the  full  cross  section  from  top-to-bot tom  and  side-to-side ,  undergoes 
no  mixing  or  longitudinal  dispersion,  and  exits  the  basin  with  a  residence 
time  equal  to  the  theoretical  residence  time.  The  theoretical  residence  time 
is  equal  to  the  ponded  volume  divided  by  the  average  flow  rate,  T  =  V/Q  . 

Such  an  ideal  regime  is  impossible  to  achieve  in  practice. 

6.  In  practice,  mixing  and  dispersion  occur  in  the  DMCA,  causing  some 


parcels  of  water  to  exit  earlier  than  the  theoretical  residence  time.  There 
are  also  "dead  zones,"  in  which  velocities  toward  the  outlet  are  considerably 
less  than  average  and  eddy  recirculation  currents  are  set  up,  causing  these 
areas  not  to  function  as  fully  involved  flow  paths  (Figure  2).  Thus,  the 
effective  flow  volume  is  less  than  the  total  DMCA  volume.  Since  the  effective 
volume  is  less  than  the  total  volume,  the  mean  residence  time  t  is  less  than 
T  . 


plan 


Figure  2.  Typical  DMCA  showing  dead  zones  and  effective  flow  volume, 

from  Montgomery  (1978) 


7.  The  mean  residence  time  of  a  DMCA  mav  be  direct  Iv  measured  by 
instantaneously  injecting  a  cjuantitv  of  tracer  dye  at  the  inflow  and  iru'niti'r- 
ing  dve  concentration  in  the  outfl('w  fFigure  1).  II  the  data  are  c(illected  in 
discrete  form,  a  set  I’f  t  ime-coticeiit  rat  ion  pairs  results.  Dye  concent  rat  ion 
dat.a  of  this  type  were  obtained  in  the  field  studies  describis!  in  this  repi'rt. 


At.  =  (t  -  t.) 

1  i+1  1 

t  h 

=  effluent  discharge  rate  during  the  i  time  interval. 

The  factors  may  be  omitted  if  flow  rates  are  constant.  Fluorescence 

readings  may  be  substituted  for  fluorescent  tracer  dye  concentration  since 
fluorescence  is  approximately  a  linear  function  of  dye  concentration.  The 
summation  ideally  should  be  continued  until  effluent  dye  concentration  returns 
to  the  background  level  and  remains  there  for  several  time  increments. 

8.  The  departure  of  residence  time  distribution  from  ideal  is  partially 
described  by  the  hydraulic  efficiency  t/T  .  The  reciprocal  of  the  hydraulic 
efficiency  is  called  the  hydraulic  efficiency  correction  factor  (HFCF) : 


HECF  =3  (2) 

t 

Required  DMCA  size  and  hence  cost  to  achieve  a  particular  effluent  quality  are 
directly  related  to  the  HFCF.  The  recommended  DMCA  design  procedure  requires 
that  the  theoretical  residence  time  as  determined  by  laboratory  settling  tests 
(or  the  DMCA  ponded  volume)  necessary  to  achieve  the  desired  effluent  quality 
be  multiplied  by  the  HECF  (Palermo,  Montgomery,  and  Poindexter  1978;  Palermo 
1985).  For  sediment  slurries  exhibiting  zone  settling,  the  design  procedure 
also  requires  that  the  minimum  DMCA  surface  area  needed  for  effective  zone 
settling  as  determined  by  laboratory  settling  tests  he  multiplied  by  the  HECF 
(Palermo,  Montgomery,  and  Poindexter  1978).  These  procedures  correct  DMCA 
sizes  determined  from  laboratorv  results  for  the  nonideal  settling  conditions 
found  in  DMCAs . 

9.  Montgomery  (1978)  recommends  assuming  a  hvdraulic  efficiency  of  0.A4 
(and  thus  a  HECF  of  2.25)  for  situations  in  which  site-specific  information  on 
the  DMCA  residence  time  distribution  is  unavailable.  Montgomery's  recommenda¬ 
tion  was  the  only  available  quantitative  design  guidance  prior  to  publication 
of  this  report.  The  impact  of  using  an  assumed  value  of  0.44  for  hvdraulic 
efficiency  when  actual  efficiencies  range  from  0.2  to  0.7  on  DMCA  sizing  cal¬ 
culations  is  shown  below; 


V  s' 


Assumed 

t/T 

Assumed 

HECF 

Actual 

t/T 

Actual 

HECF 

Pond  V'olume* 
Required 

Pond  Volume* 
from  Design 

0.44 

2.25 

0.2 

5.00 

222 

100 

0.44 

2.25 

0.3 

3.33 

148 

100 

0.44 

2.25 

0.4 

2.50 

111 

100 

0.44 

2.25 

0.44 

2.25 

100 

100 

0.44 

2.25 

0.5 

2.00 

89 

100 

0.44 

2 .25 

0.6 

1 .67 

74 

100 

0.44 

2.25 

0.7 

1.43 

64 

100 

*  Volumes  are  expressed  in  percentage  of  Che  volume  required  when 
the  actual  hydraulic  efficiency  is  0.44.  For  situations  involving 
zone  settling,  identical  numbers  would  apply  to  site  surface  area 
as  well  as  ponded  volume. 

Using  0.44  for  hydraulic  efficiency  can  significantly  increase  DMCA  size  (and 
thus  cost)  if  actual  efficiency  is  greater  than  0.44.  On  the  other  hand, 
using  0.44  might  result  in  unacceptable  effluent  quality  for  situations  in 
which  efficiencies  were  less. 

Purpose  and  Scope 

10.  The  purpose  of  this  report  is  to  provide  information  to  supplement 
the  existing  guidance  by  Montgomery  (1978)  for  selection  of  HECFs  for  DMCA 
design.  Part  11  of  this  report  describes  the  relationship  between  DMCA  flow 
patterns  and  residence  time  distribution.  A  short  literature  review  is  also 
included  in  Part  II.  Subsequent  parts  describe  collection  and  analysis  of 
residence  time  distribution  (dye  curve)  data  from  several  DMCAs.  Results  from 
DMCA  tests  are  compared  with  lab  and  prototype  tests  of  water  and  wastewater 
treatment  basins  and  are  presented  in  Part  III.  Improved  DMCA  design  guidance 
based  on  these  data  is  developed  in  Part  IV,  and  methods  for  incorporating 
this  improved  guidance  into  design  are  presented  in  Part  V. 

11.  Several  of  the  DMCA  tests  were  simulated  using  a  two-dimensional 
numerical  hydrodynamic  model.  Results  of  the  modeling  effort  are  presented  as 
Appendix  A.  Appendix  R  contains  guidance  for  collecting  DMCA  residence  time 
distribution  data,  and  Appendix  C  contains  the  DMCA  test  data. 


PART  TI:  RESIDENCE  TIME  DISTRIBUTION  AND  DMCA  HYDRAULICS 


12.  The  residence  time  distribution  of  a  DMCA  is  controlled  by  flow 
patterns  and  velocity  distributions.  This  part  contains  qualitative  descrip¬ 
tions  of  the  effects  of  hydraulic  behavior  on  residence  time  distributions  and 
the  resultant  effects  on  DMCA  suspended  solids  removal.  Relevant  literature 
is  surveyed,  including  studies  containing  measurements  of  residence  time  dis¬ 
tributions  for  DMCAs  and  similar  bodies  of  water.  Applications  of  numerical 
models  to  the  problem  of  simulating  the  hydraulics  of  DMCAs  and  predicting 
their  residence  time  distributions  are  also  briefly  reviewed. 

Departures  from  Ideal  Behavior 

13.  As  noted  previously  in  Part  I,  "ideal"  hydraulic  behavior  in  a  DMCA 
would  be  plug  flow.  Velocity  distributions  for  plug  flow  are  shown  in  Fig¬ 
ure  5.  Flow  through  real  DMCAs  strongly  departs  from  this  ideal  because  of 
unsteady  flow  rates,  wind  effects,  inlet  and  outlet  effects,  and  shear 
stresses  at  the  sides  and  bottom  of  the  DMCA.  All  of  these  effects  tend  to 
vary  continuously  in  both  time  and  space.  The  concentration  and  size  grada¬ 
tion  of  suspended  solids  in  the  Influent  dredged  material  probably  exert  an 
influence  on  hydraulic  behavior  as  well,  but  it  is  not  known  if  these  effects 
are  as  significant  as  those  produced  by  the  factors  already  mentioned.  Actual 
DMCA  conditions  are  probably  more  accurately  represented  bv  some  mix  of  advec- 
tlve  flow  and  completely  mixed  conditions,  also  shown  in  Figure  5. 

Dead  zones 

lA.  Typical  flow  patterns  for  DMCAs  are  shown  in  FiRure  6,  which  may  be 
contrasted  with  the  ideal  shown  in  Figure  3.  It  should  be  noted  that  the 
boundaries  between  the  various  zones  shown  in  Figure  b  are  vague  and  indis¬ 
tinct,  and  considerable  exchange  can  occur  between  different  "zones."  Dead 
zones  have  very  little  net  forward  flow,  although  they  may  be  either  stagnant 
or  well-mixed  Internally,  The  difference  between  an  internally  mixed  dead 
zone  and  a  well-mixed  zone  is  that  the  rate  of  exchange  of  water  between  the 
dead  zone  and  the  main  advective  flow  zone  is  much  slower. 

15.  Parcels  of  water  that  enter  dead  zones  have  verv  long  residence 
times  there,  and  a  high  percentage  of  suspended  solids  are  removed  from  them. 
However,  the  presence  of  dead  zones  adversely  affects  the  overall  treatment 


16.  Dead  zones  tend  to  occur  in  shallow  or  hydraulically  rough  areas, 
corners,  and  areas  behind  baffles  or  obstructions  that  are  sheltered  from  the 
mixing  effects  of  advectlve  flow  or  wind.  The  small  width  of  DMCA  inlets  and 
outlets  relative  to  the  width  available  for  flow  promotes  formation  of  dead 
zones  near  Inlets  and  outlets  (Figures  2,  4,  and  6).  Deposition  of  coarse 
sediments  at  the  inlet  usually  results  in  formation  of  a  fan  or  delta  as  shown 
in  Figure  3.  Such  formations  have  been  observed  at  times  to  decrease  the  size 
of  adjacent  dead  zones  by  distributing  flow  and  at  other  times  to  promote  dead 
zones  by  concentrating  flow  when  a  channel  across  the  fan  or  delta  evolves. 

17.  Dead  zones  can  also  occur  in  the  vertical  dimension,  as  well  as  the 
horizontal.  For  most  dredged  materials,  slurries  with  solids  concentrations 
greater  than  approximately  200  g/j,  do  not  readily  flow,  so  the  point  in  the 
vertical  profile  where  concentration  exceeds  200  'g/ i.  is  usually  defined  as  the 
bottom.  Slurries  with  intermediate  concentrations  (for  example,  between 

30  and  200  g/i)  flow,  but  at  distinctly  slower  rates  than  the  overlying  water, 
which  has  lower  solids  concentrations.  Therefore,  the  layer  between  the 
bottom  and  the  relatively  clarified  surface  layer  could  be  considered  a  dead 
zone  (Figure  7).  The  water  initially  contained  in  this  layer  is  trapped  for  a 
relatively  long  time.  As  the  layer  consolidates  to  concentrations  greater 
than  200  g/j,,  the  water  is  slowly  expelled  into  the  overlying  zones. 


Figure  7.  Dead  zone  in  high  concentration  layer 
near  bottom  of  DMCA 


D i spers ion 

18.  Departures  from  ideal  flow  also  occur  because  of  localized  mixing 
in  advectlve  flow,  oi  dispersion.  Dispersion  is  caused  hv  the  variation  in 
average  local  velocity.  In  general,  vel<'cities  are  higher  at  the  surface  and 
near  the  center  of  the  flow  than  at  the  bottom  or  near  the  sidi’s.  In  addi¬ 
tion,  dispersion  is  caused  bv  turbulent  mixing,  although  tliis  effect  is  lu't  as 
great  as  that  induced  bv  boundaries.  Turbulent  mixing  occurs  as  numerous. 


t  imi'-va  rv  i  ng  eddies  mix  water  in  a  given  p.arce!  with  adjacent  parcels. 


As  shown  In  the  schematic  of  ideal  flow  (Figure  5),  mixing  occurs  between 
"plugs." 

Wind  effects 

19.  Winds  can  significantly  affect  both  short-circuiting  and  dispersion 
in  DMCAs .  Alterations  in  wind  speed  and  direction  can  rapidly  modify  the  flow 
patterns  in  a  DMCA,  changing  zones  from  one  type  to  another.  Wind-induced 
surface  velocities  can  be  100  times  greater  than  mean  advective  flow  veloci¬ 
ties.  Wind-induced  surface  flows  are  accompanied  by  flows  in  the  opposite 
direction  in  lower  layers  and  spiraling  flows  in  corners  (l.iggett  and  Had- 
jitheodorou  1969).  Theoretical  calculations  by  Liu  and  Perez  (1971)  for  very 
shallow,  nonstratif led  small  water  bodies  (in  which  the  Coriolis  acceleration 
can  be  neglected)  showed  a  maximum  return  flow  velocity  of  approximately 

50  percent  of  the  maximum  surface  velocity.  The  surface  drift  direction  was 
found  to  essentially  parallel  the  wind  direction,  except  near  the  banks,  and 
to  vary  from  about  0.6  ft/sec  for  a  wind  speed  of  10  ft/sec  (6.3  mphl  to  about 
1.8  ft/sec  for  a  wind  speed  of  30  ft/sec  (20  mph) .  The  results  were  sensitive 
to  the  assumed  value  of  eddy  viscosity  (0.5  cm^/sec) ,  but  not  to  depth. 

20.  Wind  can  cause  extreme  short-circuiting  by  driving  the  entering 
flow  directly  from  the  inlet  to  the  outlet  in  less  than  5  percent  of  the  theo¬ 
retical  residence  time.  However,  observation  of  DMCAs  indicates  that  the 
maior  effect  of  wind  is  to  promote  mixing  and  not  short-circuiting.  The  high 
wind-induced  surface  velocities  and  associated  underflows  promote  lateral  and 
vertical  mixing  at  the  expense  of  dead  zones  and  advective  flow  zones.  Thus, 
wind  tends  to  increase  the  fraction  of  the  DMCA  dominated  by  compl  e  t  e  1  y-m.ixed 
cond i t ions . 

Dye  Tracer  Curves  and  Residence  Time  Distributions 

Mean  residence  time 

21.  As  noted  previously  in  Part  I,  the  curve  of  DMCA  effl\)ent  dve  c('n- 
centrntion  versus  time  that  results  from  a  slug  input  of  dve  into  the  inflow 
is  identical  to  a  frequency  distribution  of  residence  times.  The  presence  ol 
dead  zones  is  manifest  in  the  dve  tracer  curve,  the  mean  residence  tine  t 
being  less  than  tlie  theoretical  residence  time  T  .  The  mean  resideiue  time 
mav  be  computed  from  the  dve  tracer  curve  \ising  Fqiiatii'ii  1.  The  li  imens  icni  1  ess 


1  ^ 


ratios,  hydraulic  efficiency  t/T  ,  and  HK(!F  T/t  are  numerical  measures  of 
Che  influence  of  dead  zones. 

Dispersion  index 

22.  Gross  longitudinal  mixing  is  manifested  in  the  dve  tracer  curve  bv 
the  spread  of  the  base  of  the  curve.  The  dispersion  index  d  is  a  dimension¬ 
less  number  that  measures  the  amount  of  mixing  occurring.  The  dispersion 
index  is  obtained  by  dividing  the  variance  of  the  residence  time  distribution 
by  the  mean  residence  time  squared. 


where 


=  variance  of  the  residence  rime  distribution 


t "  =  mean  residence  time  squared 

The  dispersion  index  mav  be  computed  t rom  a  digital  dve  tracer  :;urve  using  the 
following  formula,  where  previous  definitions  of  vari.ililes  .ipplv: 
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A  high  value  for  tlie  d  i  spe  r  5;  i  (tn  index  indicates  a  wide  spreati  ol  i  I  ow- 1  h  rt'u  ih 
times.  Also,  because  of  continuitv.  a  high  dispersiiui  index  indiiate<  that  a 
large  tr.iction  of  t  lie  flow  exit-^  earlier  than  t  . 

hfficiencv,  disper¬ 
sion,  and  sol  ills  removal 

2  5.  Ihe  dispersion  imiex  and  the  mean  re'-ideiue  time  ari-  entire!'.' 


independent  of  one  another.  For  example,  l  igiita-  ■''.i  sh. 
identical  dispersion  indices  hut  ditterent  me, in  it'-^i.ie: 
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solids  removal  efficiency  because  of  the  nonlinearity  of  the  solids  removal 
curve  (Figure  1).  The  greater  percentage  of  removal  of  solids  from  parcels  of 
water  with  residence  times  greater  than  t  is  not  great  enough  to  counterbal¬ 
ance  the  lesser  percentage  of  removal  of  solids  from  parcels  of  water  with 
residence  times  less  than  t  .  Figure  9  contrasts  the  characteristics  of  two 
extreme  types  of  residence  time  distributions:  plug  flow  and  complete  mix. 


PLUG  FLOW 

t/T  =  1 
d  =  0 


COMPLETELY  MIXED 

t/T=  1 
d  =  1 


t/T 

Figure  9.  Residence  time  distributions  from  plug  flow 
and  completely  mixed  basins 


Conceptual  model 

24.  A  simple  conceptual  model  may  be  used  to  reproduce  residence  time 
distributions  similar  to  that  produced  by  the  complex  situations  shown  in 
Figure  6.  Figure  10  shows  a  conceptual  model  composed  of  four  compartments: 
an  advective  zone  at  the  inlet,  and  then  a  completely  mixed  zone,  an  advective 
zone,  and  a  dead  zone  to  represent  the  main  body  of  the  DMCA.  As  can  be  seen 
from  Figure  10,  the  overall  DMCA  residence  time  distribution  is  the  sum  of  the 
outputs  from  the  three  main  compartments  and  approximates  a  log-normal  distri¬ 
bution.  The  exact  shape  of  the  log-normal  distribution  depends  on  the  disper- 
sioi.  of  the  output  from  the  four  compartments  and  the  relative  magnitude  of 
the  compartment  voltimes  and  flow  rates.  Figure  10  is,  of  course,  a  simplified 
analogy.  Some  DMCAs  mav  have  small  complete-mix  zones  in  l.irge  scour  holes  at 
the  point  of  inflow,  instead  of  advective  flow  fAj).  Cons iderable  exchange 
that  is  not  part  of  the  Figure  10  .an.ilogv  occurs  hetweei'  well-mixed, 
advective,  and  dead  zones. 

2  S .  The  rel.ative  volumes  of  the  f('iir  c(imp,artment  s  .and  tlie  division  I'f 
flow  amont’  the  1, at  ter  three  v.arv  c  mi  t  i  nuous  1  v  ,as  functions  of  w  i  ml  ,  ti'tal 
flow,  DMCA  I'otfom  t  opi’g  r  .aph  v ,  ,an<!  pl.an  geometr\',  eti'.,  so  quant  i  t  .a  t  i  ve  us,,  of 
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Reported  Residence  Time  Distributions 


28.  The  literature  contains  little  Information  about  the  hydraulics  of 
large,  shallow  ponds,  such  as  DMCAs,  relative  to  the  volume  of  Information 
about  sedimentation  basins  for  water  and  wastewater  treatment.  This  situation 
is  partially  attributable  to  the  difficulties  of  performing  tracer  tests  on 
large  facilities  with  long  residence  times  and  remote  locations.  Relevant 
residence  time  distribution  studies  are  summarized  below. 

DMCA  studies 

29.  Montgomery  (1978)  reported  the  results  of  three  dye  tests  at  two 
DMCAs.  In  all  three  tests,  short-circuiting  was  significant,  and  the  hydrau¬ 
lic  efficiency  for  all  three  tests  was  reported  as  less  than  0.5. 

30.  Poindexter  and  Perrier  (1980)  reported  the  results  of  eight  DMCA 
dye  tests,  including  the  three  previously  reported  by  Montgomery  (1978). 
Results  of  the  dye  tests  were  sensitive  to  the  unsteady  nature  of  dredge  dis¬ 
charge  and  wind.  The  sites  studied  had  low  values  of  length-to-width  ratio 
(L/W)  (2  to  4) ,  and  none  of  the  numerical  descriptors  of  the  dye  curves  were 
found  to  be  correlated  with  L/W  ratio,  however,  errors  were  made  in  dye  curve 
analysis,  leading  to  slightly  incorrect  values  for  the  statistical  descrip¬ 
tors.  Wind  caused  serious  short-circuiting,  particularly  in  shallow  DMCAs. 
WTiile  wind  affected  the  shapes  of  the  dye  curves,  mean  residence  times  were 
less  affected. 

31.  Palermo  (1984)  reported  the  results  of  three  dye  tests  conducted  in 
conjunction  with  his  research  on  pollutant  removal  in  dredged  material  con¬ 
tainment  areas.  The  primary  purpose  of  the  tests  was  to  estimate  the  mean 
residence  time  so  that  effluent  sampling  could  be  lagged  from  influent  sam¬ 
pling  by  approximately  that  time.  However,  only  two  of  the  tests  (at  Mobile 
Harbor  and  Norfolk  Harbor")  produced  complete  tracer  curves.  estimates  of 
average  flow  rate  and  DMCA  dimensions  were  not  adequate  for  accurate  determi¬ 
nation  of  T  ,  and  thus  t/T  . 

32.  Gallagher  and  Company  (1978)  produced  a  report  specificallv  on  DMCA 
residence  time  distributions.  Thev  presented  results  of  four  dve  tests,  one 
of  which  was  obtained  from  Montgomery  (1978).  On  the  basis  of  these  tests, 
they  concluded  that  "...  considerable  short-circuiting  occurs  in  large  shallow 
basins,  and  this  results  in  inefficient  flow  patterns.  Actual  retention  times 
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are  significantly  shorter  than  ideal  through-flow  times,  and  wind  effects  can 

be  very  detrimental," 

Water  and  wastewater 
treatment  basin  studies 

33.  Marske  and  Boyle  (1973)  reported  the  results  of  dye  tests  on  proto¬ 
type  chlorine  contact  chambers.  They  tested  a  wide  variety  of  basin  shapes 
and  baffling  arrangements.  Theoretical  residence  times  for  the  rectangular 
basins  ranged  from  0.46  to  7.78  hr.  Flow  rates  ranged  from  0.46  to  47  cfs. 

34.  Marske  and  Boyle  (1973)  ran  a  pair  of  tests  with  identical  hydrau¬ 
lic  conditions,  but  with  a  20-mph  wind  blowing  directly  toward  the  outlet  in 
one  case  and  directly  away  from  the  outlet  in  the  other  case.  Wind  was  found 
to  cause  serious  short-circuiting.  The  ratio  of  initial  time  of  dye  appear¬ 
ance  to  T  decreased  from  0.54  with  the  wind  blowing  upstream  to  0.36  with 
the  wind  blowing  downstream.  However,  mean  residences  for  the  two  tests  were 
essentially  equal. 

35.  The  dispersion  index  d  was  recommended  as  the  most  reliable  indi¬ 
cator  of  basin  hydraulic  performance.  A  value  of  L/V  greater  than  20  was 
necessary  to  reduce  the  dispersion  index  below  0. 1 ,  an.,  a  value  of  L/W  greater 
than  40  was  necessary  to  achieve  the  minimum  value  of  d  ,  which  was  about 
0.02.  Longitudinal  baffling  was  found  to  be  better  than  transverse  baffling, 
because  a  larger  value  of  L/W  can  be  achieved  with  fewer  baffles.  In  addi¬ 
tion,  fewer  baffles  result  in  fewer  corners  or  shadow  areas,  which  are  poten¬ 
tial  dead  spaces, 

36.  Mangelson  and  Watters  (1972)  conducted  tracer  tests  on  prototype 
waste-treatment  lagoons  at  Logan,  I'tah,  and  in  the  laboratory  on  a  20-  by  40- 
by  3.5-ft  model.  Various  inlet  and  outlet  conditions  were  used  in  the  model 
tests.  The  average  hydraulic  efficiency  for  the  lagoons  was  about  0.58.  The 
lagoons  had  values  of  L/W  between  l.b  and  2.3.  Values  of  hvdraulic  efliciency 
for  the  unbaffled  model  ranged  from  0.52  to  0.64  (omitting  the  highest  and 
lowest  values)  for  depths  of  1.5  ft  and  L/W  equal  to  2.0.  In  model  tests  with 
baffles  in  various  arrangements,  a  value  of  I./W  of  5.26  produced  a  hvdraulic 
efficiency  of  0.72,  and  values  ot  L/W  from  13  to  45  produced  values  (if  hvdrau¬ 
lic  efficiency  between  O.Hl  and  O.Hf),  except  for  one  experiment  in  which  the 
battles  did  not  overlap.  In  this  case,  the  value  ol  lo.'draulic  etliciencv  was 


37.  In  all  of  the  prototype  tests,  dimensionless  initial  time  t  . /T 

was  less  than  0.1,  and  dimensionless  time  of  peak  concentration  t  /T  was 

P 

between  0.1  and  0.2,  denoting  extreme  dispersion.  For  the  unbaffled  model 

tests,  t./T  was  between  0.06  and  0.18,  and  t  /T  was  between  0.18  and  0.30. 

^  P 

The  introduction  of  baffles  decreased  the  amount  of  dead  space  and  signifi¬ 
cantly  increased  values  of  t  /T  and  t./T  .  In  general,  the  more  baffles 

pi 

there  were  and  the  closer  the  baffles  came  to  the  opposite  wall,  the  closer 
the  residence  time  distribution  approached  plug  flow. 

Mathematical  Models 

38.  Gallagher  and  Company  (1978)  developed  a  potential  flow  model 
(laminar  flow)  and  applied  it  to  hypothetical,  rectangular  DMCAs  with  various 
inlet  and  outlet  configurations  and  spur  dike  arrangements.  Although  the 
model  is  not  a  realistic  depiction  of  DMCA  conditions,  cannot  be  considered 
quantitative,  cannot  show  flow  separation  or  dead  zones  (although  areas  of 
very  low  predicted  velocity  can  be  interpreted  as  dead  zones)  ,  and  cannot 
handle  wind  effects,  it  is  a  qualitative  indicator  of  hydraulic  behavior. 

Based  on  model  output,  Gallagher  and  Company  recommended  outlet  weirs  as  long 
as  possible  and  one  to  four  baffles  extending  75  percent  of  the  way  across  or 
along  the  basin.  Economics  of  spur  dikes  were  discussed,  and  tables  were 
presented  showing  the  increase  in  L/W  ratio  and  dike  construction  cost 
associated  with  various  spur  dike  arrangements. 

39.  Koussis,  Saenz,  and  Thackston  (1982)  evaluated  available  hydrody¬ 
namic  models  with  respect  to  their  ability  to  reproduce  the  important  phe¬ 
nomena  occurring  in  DMCAs.  Fconomic  considerations  associated  with  the  use  of 
models  for  design  were  also  evaluated.  Koussis,  Saenz,  and  Thackston  sug¬ 
gested  development  of  a  steady-state,  two-dimensional,  vertically  integrated 
hydrodynamic  model  capable  of  handling  natural  topographies,  realistic  bound¬ 
ary  conditions,  and  moderate  wind  effects.  Such  a  model  could  be  coupled  with 
a  sediment  transport  model  for  prediction  of  DMCA  effluent  solids  concentra¬ 
tion.  Thev  recognized  that  the  two-dimensional  hydrodynamic  model  might  fail 
to  adequately  describe  DMCA  hydrodynamics,  and  in  that  case  recommended 
development  of  a  three-dimensional  model. 

40.  In  Appendix  A,  Bach  reports  the  results  of  an  effort  to  simulate 
the  hydraulic  behavior  of  DMCAs  using  three-dimensional  and  two-dimensional 


computer  models.  Satisfactory  results  were  obtained  from  initial  runs  of  the 
three-dimensional  model,  but  the  modeling  effort  was  continued  with  the  two- 
dimensional  codes  to  meet  the  time  and  money  constraints  of  this  study.  The 
two-dimensional  model  was  vertically  averaged,  and  thus  it  could  not  simulate 
circulation  in  the  vertical  plane  caused  by  wind  shear.  Nevertheless,  the 
model  was  able  to  predict  peaks  and  means  of  eight  observed  DMCA  residence 
time  distributions  when  model  coefficients  were  chosen  bv  trial  and  error. 

41.  Selection  of  model  coefficients  for  predictive  use  of  the  model 
remains  problematic,  even  though  a  linear  regression  equation  for  the  model 
coefficients  based  on  surface  drift  velocity,  DMCA  surface  area,  discharge, 
and  wind  direction  was  derived.  Since  only  eight  data  points  were  used  in 
fitting  the  equation  (which  has  four  independent  variables)  and  since  tlie 
field  data  used  to  adjust  the  model  are  of  low  quality,  use  of  the  regression 
equation  for  selecting  model  coefficients  for  other  DMCAs  is  not  generally 
recommended . 

42.  Tatom  and  Waldrop  (1986)  reported  Initial  work  on  the  adaptation  of 
a  vertically  averaged,  two-dimensional  computer  model  for  simulation  of  ash 
settling  ponds  associated  with  coal-fired  power  plants.  Pond  geometries  are 
similar  to  large  DMCAs,  and  flow  rates  are  slightly  lower.  The  model  will 
Include  inputs  of  water  to  the  pond  from  rainfall  and  surface  runoff.  Plans 
for  calibration  include  collection  of  extensive  field  data.  Results  of  this 
work  may  be  useful  to  future  DMCA  designers. 

Summary  of  Literature 

43.  Little  information  is  available  regarding  residence  time  distribu¬ 
tions  of  large,  shallow  impoundments.  Field  data  are  scarce  because  of  the 
extremely  long  residence  times  and  attendant  time  requirements  for  dve  tests. 
In  general,  these  impoundments  tend  to  have  ef f ic ieiic  ies  in  the  0.4  to 

0.6  range  and  high  dispersion  indices.  W'ind-induced  currents  can  be  two 
orders  of  magnitude  greater  than  average  advective  currents. 

44.  Callagher  and  Company  (1978)  present  qualitative  guidance  for 
hydraulic  design  of  DM(L\s.  Thev  discuss  the  benefits  and  costs  of  spur  dikes. 
Test  data  from  chlorine  contact  chambers  and  physical  models  sluiw  that 


dispersion  and  short-circuiting  may  be  greatly  reduced  by  increasing  L/W 
(which  may  be  done  with  spur  dikes  or  baffles*). 

45.  Generation  of  residence  time  distributions  by  simulating  DMCA 
hydraulics  using  computer  models  appears  to  be  just  within  the  state  of  the 
art.  However,  the  costs  of  collecting  adequate  field  data  for  calibration  and 
running  the  models  are  obstacles  to  their  use.  Vertically  averaged,  two- 
dimensional  models  are  not  able  to  simulate  wind-induced  circulation  in  the 
vertical  plane.  Field  studies  indicate  that  wind-induced  circulation  affects 
the  overall  shape  of  the  residence  time  distribution  (reflected  in  parameters 
like  the  time  to  peak  and  the  dispersion  index)  but  not  the  mean  residence 
time  or  hydraulic  efficiency. 


*  The  term  "spur  dikes"  ns  used  herein  refers  to  baffles  in  a  DMC^A .  DMCAs 
are  sometimes  baffled  hv  ccinstruct ing  interior  earthen  embankments.  Water 
and  wastewater  treatment  basins  are  usual Iv  baffled  with  reinforced  concrete 
walls.  Wooden  baffles  mav  be  used  in  moilel  tests.  The  use  of  floating 
baffles  in  DMCAs  is  discussed  in  Fart  V. 


PART  III:  DESCRIPTION  OF  FIELD  TESTS 


Data  Sources 


46.  Records  of  13  DMCA  dye  tracer  tests  were  available.  All  of  the  dve 
tracer  tests  were  conducted  as  part  of  the  US  Army  Engineer  Waterways 
Experiment  Station  (WES)  research  projects  at  active  CE  dredging  sites.  In 
most  cases,  the  tests  were  intended  to  provide  data  to  support  research  on 
subjects  other  than  residence  time  distribution  or  hydraulic  efficiency. 
Therefore,  the  quantity  and  quality  of  data  from  the  various  dye  tests  varv 
widely.  In  an  attempt  to  standardize  the  results  for  analysis,  all  available 
original  data  were  reexamined  and  reanalyzed  in  a  consistent  manner. 

47.  One  of  the  13  tests  had  to  be  eliminated  because  the  dye  tracer 
curve  was  incomplete  and  because  some  essential  physical  data  were  missing. 

The  other  12  tests  are  listed  in  Table  1.  Normalized  dye  curves,  site  maps, 

and  wind  data  are  displayed  on  similar  scales  for  all  12  tests  in  Fig¬ 

ures  12  through  23.  Dye  concentrations  were  normalized  by  dividing  by  C  , 

3T 

C„.^/Cdt  (5) 

o 

Physical  Variables 

48.  The  physical  variables  of  length,  width,  area,  and  volume  were  all 

recalculated,  starting  with  the  map  of  the  site.  These  maps  were  digitized, 

and  the  areas  were  calculated.  Smooth  flow  paths  from  inlet  to  outlet  were 

sketched,  and  their  lengths  were  measured  and  used  as  1.  in  subsequent 
analyses.  Widths  were  measured  as  the  distance  from  boundary  to  boundary 
perpendicular  to  the  flow  path.  Several  typical  widths  were  averaged  to  a 
single  value. 

49.  The  depths  of  the  basins  were  measured  bv  vari<Mis  methods.  Mont¬ 
gomery  (1978)  took  samples  at  various  depths,  analyzed  each  tor  suspended 
solids,  plotted  the  data,  and  defined  depth  as  the  distance  f r(im  the  w.iter 
surface  to  the  point  where  the  solids  concent  rat  ion  was  about  20(1  g/'.  I’oin- 
dcxter  and  Perrier  (1980)  used  a  hydrometer  consisting  ot  a  small  plastic  bot¬ 
tle  filled  with  a  sediment-water  mixtiire,  adjusted  so  that  it  was  neutrally 
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Table  1 


hocat  ion  of 

Tests  and  Data 

Sources* 

Test 

Site 

Date 

R  e  1  e  r  e  n  c  e 

! 

Yazoo  River  '15,  Miss. 

2/23/77 

(lallagher  and  Co.  (19 /S') 
Mont  gomerv  ( I 9  78) 

n 

Fort  Fust  is,  Va . 

5/12/79 

Poindexter  and  Perrier 
(1980) 

3 

Burnsvi lie,  Miss. 

3/81 

Schroeder  (1983  "I 

4 

Yazoo  River  f'S,  Miss. 

3/17/77 

Montgomerv  (1978’) 

5 

Y.azoo  River  fl5.  Miss. 

8/9/78 

Poindexter  and  Perrier 
(19801 

6 

Mobile  (N.  Blakelyl,  Ala. 

6/29/78 

Poindexter  and  Perrier 
(1980) 

7 

Yazoo  River  /16,  Miss. 

2  or  3/77 

Callagher  and  Co.  (IDfR) 
Poindexter  and  Perrier 
(19801 

8 

Yazoo  River  #3,  Miss. 

12/12/78 

Poindexter  and  Perrier 
(19801 

9 

Mobile  (N.  Blakely),  Ala. 

7/7/82 

Palermo  (1984'i 

11 

Norfolk  (Craney  Is.),  Va. 

2/7/83 

Pa  1  (>rmo  (  I  984 1 

12 

Black  Rock  Harbor  (I'nited 
Illuminating),  Conn. 

10/26/83 

I'npublished  notes,  site 
description  bv  P.ilermo 
and  Pr.inger  (19841 

13 

Fowl  River,  Ala. 

4/i:’/77 

Montgcxmerv  (19^81 

*  Test  10  was  not  used  because  the  dve  curve  was  incomplete  ami  s  i  1 1’  peometrv 
data  were  poor. 


buoyant  when  immersed  in  a  slurrv  with  a  tc>tal  solid;;  (npcentr.it  ion  ni 
200  g/ P  .  This  was  I(/wered  into  the  basin  bv  a  ci'rd  (  rom  a  I'oat  until  it 
stopped,  and  the  depth  was  recorded  as  the  dist.ance  from  the  center  ot  the 
bottle  to  the  water  surface.  Schroeder  fl'tH'n  estimated  depths  h.ised  op  t  i|e 
difference  fietween  me.i.sured  water  surface  elev.it  ion  and  elevat  iom  o'  the 
bottom  taken  t  rom  as-built  plans.  Schroedet's  test  occurted  e.itiv  t!.e  dis- 


p('’s,il  operation  before  much  deposition  had  I'ccut  t  cd  in  the  h,i‘  in. 


WINDDIPtCTION  DEGRtES 


P'igure  13.  Test  2,  Fcrt  Fust  is,  Virgini.i.  (al  W  i  lu!  d  i  rec  t  i  ui’ .  Y-axis  sln'ws 
liirection  from  which  the  wind  blew  in  <legrees  clockwise  1  rc'u-  north,  i.e.,  'hV’ 
for  east,  180°  for  south,  270°  for  west.  A  wind  direct  ioic  oi  O'’ 

!T 


indicates  calm.  f  b  I  Wind  speed.  (c)  Norma  I  i  j-ed  dve  curve.  ( 


1 


f  r  dt . 


Missing  valuer:  of  concent  rat  ion  at  large  t  ''I'  estimated  using  pis'cc'dure 
described  in  I’art  V.  cl)  Site  plan,  1  in.  =  1,000  ft.  Matcliing  indicates 

drv  area 


cmdS  QNI.V 


!  T 


?'igure  15.  Test  4,  Yazoo  River  No.  5,  Mississippi.  (a)  Wind  direction. 
Y-axis  shows  direction  from  which  the  wind  blew  in  degrees  clockwise  from 
north,  i.e.,  90°  for  east,  180°  for  south,  ?70°  for  west.  A  wind  direc¬ 
tion  of  0°  indicates  calm.  (h)  Wind  speed.  (c)  Normalized  dve  curve. 

I 

“  y  y*  C  •  Missing  values  of  concentration  at  large  t/T  estimated 
o 

using  procedure  described  in  Part  V.  (d)  Site  plan,  1  in.  =  1,000  ft. 

Hatching  Indicates  drv  area 


Figure  16.  Test  ,  Ynzno  River  No.  6,  Mississippi.  <'a)  Wind  direction. 
Y-axis  shows  direction  from  which  the  wind  blew  in  decrees  clockwise  from 
n(’rth,  i.e.,  60"  for  east,  IHO”  for  s(nith,  270"  for  west.  A  wind  direc¬ 
tion  ot  0^  indicates  caliTi.  (hi  Wind  speed.  (cl  Normalized  dve  curve. 

1 

(■  =  r  i'  (it.  Missing  Veilues  cf  concent  r.U  i  on  nt  large  t /T  estimated 

I 

(1 

usin^;  procedure  described  in  Part  V.  ( d  1  Site  plan,  I  in.  = 

Hatching  indicates  drv  area 


1 ,000  ft  . 


Figure  17.  Test  6,  Mobile  (N.  Blakely),  Alabama.  (a)  Wind  direction. 
Y-axis  shows  direction  from  which  the  wind  blew  in  degrees  clockwise  from 
north,  i.e.,  90”  for  east,  180”  for  south,  270”  for  west.  A  wind  direc¬ 
tion  of  0”  Indicates  calm.  (b)  Wind  speed.  (cl  Normalized  dve  curve. 

1 

C  =  y  y"  0  dt.  Missing  values  of  concentration  at  large  t/T  estimated 
o 

using  procedure  described  in  Part  V.  (d)  Site  plan,  1  in.  =  1,000  ft. 

Hatching  indicates  drv  area 


•■w-V 


DATA  NOT 
AVAILABLE 


DATA  NOT 
AVAILABLE 


Figure  18.  Test  7,  Yazoo  River  No.  6,  Mississippi.  (a)  Wind  direction. 
Y-axis  shows  direction  from  which  the  wind  blew  in  degrees  clockwise  from 
north,  i.e.,  90°  for  east,  180°  for  south,  270°  for  west.  A  wind  direc¬ 
tion  of  0°  indicates  calm.  (b)  Wind  speed.  (c)  Normalized  dve  curve. 

1 

=  —  y  C  dt.  Missing  values  of  concentration  at  large  t/T  estimated 
o 

using  procedure  described  in  Part  V.  (d)  Site  plan,  I  in.  =  1,000  ft. 

Hatching  indicates  drv  area 


Figure  19.  Test  8,  Yazoo  River  No.  3,  Mississippi.  (a)  Wind  direction. 
Y“axis  shows  direction  from  which  the  wind  blew  in  degrees  clockwise  f rom 
north,  i.e.,  90°  for  east,  180°  for  south,  270°  for  west.  A  wind  direc¬ 
tion  of  0  indicates  calm.  (b)  Wind  speed.  (cl  Normalized  dve  curve. 

1  V 

~  y  0  dt.  Missing  values  of  concentration  at  large  t/T  estimated 
o 

using  procedure  described  in  Pari  V.  (d)  Site  plan,  I  in.  =  1,000  ft. 

Hatching  indicates  dry  area 


Figure  20.  Test  9,  Mobile  (North  Blakely),  Alabama.  (a)  Wind  direction. 
Y-axis  shows  direction  from  which  the  wind  blew  in  degrees  clockwise  from 
north,  i.e.,  90°  for  east,  180°  for  south,  270°  for  west.  A  wind  direc¬ 
tion  of  0°  indicates  calm.  (b)  Wind  speed.  (c)  Normalized  dve  curve. 

3T 


C 


1 


y*  C  dt.  Missing  values  of  concentration  at  large  t/T  estimated 


using  procedure  described  in  Part  V.  (d)  Site  plan, 

Hatching  indicates  drv  area 


1  in.  =  1 , 000  f  t , 


WIND  SPEED,  MPH 


;’OUO  'I  .’(K)o  i  I 


Figure  21.  Test  II,  Norfolk  (Craney  Island),  Virginia.  (a)  Wind  direction. 
Y-axis  shows  direction  fron  which  the  wind  blew  in  degrees  clockwise 
from  north,  i.e.,  90°  for  east,  180°  for  south,  270°  for  west.  A  wind 
direction  of  0°  indicates  calm.  (b)  Wind  speed.  (c)  Normalized  dye  curve. 

1 

C  =  —  r  C  dt.  Missing  values  of  concentration  at  large  t/T  estimated 

o  1  •' 

o 

using  procedure  aescribed  in  Part  \' .  (d)  Site  plan,  I  in.  =  6,666  ft. 

Hatching  indicates  dry  area 


■  ,•••  w*-  .*•  .*•  - 


Urn 


soo  0  600  rr 


Figure  22.  Test  12,  Black  Rock  Harbor  (United  Illuminating),  Connecticut, 
(a)  Wind  direction.  Y-axis  shows  direction  from  which  the  wind  blew  in 
degrees  clockwise  from  north,  i.e.,  90°  for  east,  180°  for  south,  270°  for 
west.  A  wind  direction  of  0°  indicates  calm.  (b)  Wind  speed. 

(c)  Normalized  dry  curve. 


1 

C  =  —  y*  C  dt.  Missing  values  of  concentration  at  large  t/T 
o 


estimated  using  procedure  described  in  Part  V.  (d)  Site  plan,  1  in.  = 
1,000  ft.  Hatching  indicates  dry  area 


Figure  23.  Test  13,  Fowl  River,  Alabama.  (a)  Wind  direction.  Y-axis  shows 
Y-axis  shows  direction  from  which  the  wind  blew  in  degrees  clockwise 
from  north,  i.e.,  90°  for  east,  180°  for  south,  270°  for  west.  A  wind 
direction  of  0°  indicates  calm.  (b)  Wind  speed.  (c)  Normalized  dye  curve, 
j  3T 

C  =  —  y  C  dt.  Missing  values  of  concentration  at  large  t/T  estimated 
o 

using  procedure  described  in  Part  V.  (d)  Site  plan,  1  in.  =  1,000  ft. 

Hatching  indicates  dry  area 


50.  Volumes  were  calculated  by  multiplying  area  times  average  depth. 

The  recalculated  physical  variables  and  dimensionless  ratios  for  each  test  are 
shown  in  Table  2.  Most  of  the  sites  were  approximiitely  rectangular,  with  I,/W 
ratios  varying  from  0.7  to  4.1.  Most  h/W  ratios  were  between  2.0  and  3.0. 

All  sites  except  the  one  at  Fort  Fust  is  were  unbaffled. 


Table  2 

Phvsical  Variables  at  Fach  Test  Site 


Test 

Length 

Width 

Depth  Area 

Vol ume 

No. 

ft 

ft 

f  t  acres 

acre- f  t 

L/W 

W/D 

I./D 

1 

I  ,680 

410 

7.4  15,4 

1  lO.O 

4  .  1 

5  5 

2  30 

O 

1,550 

570 

2.7  23.4 

63.0 

2.7 

210 

5  70 

3 

450 

140 

2.7  1.7 

4.5 

3.2 

52 

I  70 

4 

1  ,090 

410 

2.0  10.3 

21.0 

2.7 

210 

540 

5 

2,000 

610 

8.0  28.4 

2  30.0 

3.3 

76 

250 

6 

3,010 

1  ,350 

6.4  91,5 

586.0 

O  O 

210 

470 

7 

1,120 

460 

6.0  11.7 

70.0 

2.4 

77 

190 

8 

2,540 

750 

7.0  34.3 

240,0 

3.4 

107 

360 

9 

1.L 

12 

1,850 

1,320 

56.3 

- 

1  .4 

- 

- 

220 

140 

3,2  0.7 

2.3 

1.6 

44 

68 

13 

820 

410 

2.2  10,7 

2.4 

2.0 

190 

3  70 

51.  The 

inlet  device 

in  all  cases  except  Test  3  was 

a  singl 

e  pipe 

0  f 

the  same  size 

as  the  dredge 

discharge  pipe. 

Test  3  was  conducted 

in  the 

secondary  cell 

of  a  two-celled  disposal  area 

.  The  inlet  was  the  corrugated 

metal 

culvert 

that  released 

water  from  the  box-type  weir  draining 

the  primarv 

area . 

1 1  was 

the  only  submerged  inlet. 

52.  In 

most  cases,  the  overflow  weir 

was  relativelv 

short  c 

ompared 

with 

the  basin  width,  and  in  all 

cases  except  one 

( Fowl  River) , 

i  t  wa  s 

at  the 

oppo- 

site 

end  of  the  long  dimension  of  the  basin 

f rom  the  inlet 

pipe. 

In  eight 

tests 

,  the  outlet  device  was  a  straiglit  weir 

100  to  450  ft 

in  length.  In 

three  cases,  a  box-type  drop  inlet  was  used. 


Flow  Rates 


53.  Inflow  rates  were  inconsistently  estimated.  In  no  case  was  the 
actual  pipeline  flow  rate  measured.  Schroeder  (1983)  estimated  the  outflow 
rate  for  Test  3  by  measuring  the  head  on  the  overflow  weir.  Montgomery  ('1978) 
determined  inflow  by  multiplying  the  cross-sectional  area  of  the  inflow  pipe 
by  an  estimated  velocity,  ranging  from  12  to  20  ft/sec,  and  checked  these 
figures  by  estimating  outflow  by  measuring  the  head  on  the  overflow  weir 
(Tests  I  and  4).  In  all  three  cases,  the  flow  was  essentially  constant  and 
uninterrupted  during  tests. 

54.  Poindexter  and  Perrier  (1980)  did  not  report  the  basis  for  their 
flow  estimates.  The  estimated  mean  pipeline  velocities,  calculated  from 
reported  flow  rates  and  pipe  sizes,  varied  from  10.5  to  18.5  ft/sec  for  test 
durations  of  80  to  339  hr.  For  a  dredge  to  operate  uninterrupted  for  this 
long  is  inconceivable.  Montgomery  (1978)  analyzed  operation  records  for  two 
dredges  and  found  that  average  operation  times  were  15  to  17  hr/day.  There¬ 
fore,  Poindexter  and  Perrier's  flow  rates  were  recomputed  based  on  pipeline 
sizes,  an  average  velocity  of  15  ft/sec,  and  16  hr/day  of  operation. 

55.  The  flow  rate  at  the  North  Blakely  site  (Test  9)  was  found  by  using 
records  of  dredge  operation  to  determine  the  percentage  of  the  dye  test  dura¬ 
tion  that  the  dredge  was  pumping  and  by  assuming  a  pipeline  velocity  of 

12  ft/sec.  This  lower  value  of  velocity  is  reasonable  in  light  of  the  facts 
that  the  pipeline  was  about  13,000  ft  long  and  there  was  no  booster  pump.  The 
flow  rate  at  the  Black  Rock  Harbor  test  was  estimated  by  measuring  the  head  on 
the  overflow  weir.  Gallagher  and  Company  (1978)  did  not  report  the  flow  rate 
at  their  test  at  the  Yazoo  River  No.  6  site,  but  Poindexter  and  Perrier  (1980) 

3 

reported  it  as  20.4  ft  /sec.  This  was  assumed  to  be  correct. 

56.  The  inlet  pipes  varied  in  diameter  from  4  in.  (Black  Rock)  to 

27  in.  (North  Blakely),  but  five  were  18  in,,  and  two  were  16  in.  All  dredges 
were  floating  cutterhead  dredges,  except  at  Black  Rock  Harbor.  There,  a  clam¬ 
shell  dredge  filled  two  barges,  and  material  was  pumped  from  the  barge  by  a 
submersible  pump  through  a  4-in.  flexible  hose.  The  pipe  diameters,  flow 
rates,  and  types  of  dredged  material  are  tabulated  in  Table  3. 


Pipe  Dlam 
In . 


Flow  rate 
f  t^/sec 


Theoretical 
Residence 
Time ,  V/Q, *  hr 


Duration** 

hr 


Type  of 
Dredged 
Material 


respectively.  Wind  speeds  exhibited  a  modest  positive  correlation  (r 

2 

=  0.252),  but  wind  directions  did  not  (r*"  =  0.028).  Scatter  plots  showed 
definite  associations  of  both  speed  and  direction  at  the  two  stations,  but 
correlations  were  influenced  by  extreme  outliers.  The  Greenwood  measurements 
were  used  to  calculate  the  average  wlndspeeds  used  in  the  analysis  below. 

Also,  wind  data  in  Figures  12,  15,  16,  and  19  are  from  Greenwood.  Slade 
(1968)  noted  that  wind  measurements  at  adjacent  stations  in  flat  terrain  show 
high  levels  of  association,  but  spatial  variation  is  much  greater  in  rugged 
terrain. 

59.  Analysis  of  wind  speed  and  direction  data  was  hampered  by  a  lack  of 
data  regarding  test  starting  times.  No  dates  were  found  for  Test  7,  and  dates 
but  not  times  were  found  for  Tests  I,  4,  5,  6,  and  7.  A  summary  of  wind  data 
is  in  Table  4.  Most  of  the  wind  statistics  in  Table  4  are  based  on  observa¬ 
tions  at  3-hr  intervals.  When  exact  test  starting  times  were  unavailable, 
reasonable  times  were  assumed.  Resultant  wind  directions  (the  direction  from 
which  the  wind  blows)  were  determined  by  summing  the  wind  vectors  measured  at 
3-hr  Intervals.  Observations  at  3-hr  intervals  were  not  available  for  any 
locations  in  close  proximity  to  Test  3. 

60.  Surface  drift  velocities  were  computed  at  3-hr  intervals  using  the 
method  described  by  Gallagher  and  Company  (1978).  Gallagher  and  Company 
(1978)  gave  equations  for  surface  drift  velocity  as  a  function  of  wind  speed 
and  eddy  viscosity.  A  power  function  was  fitted  to  Kaurashev's  data  (1960,  in 
Gallagher  and  Company  1978)  to  express  eddy  viscosity  as  a  function  of  depth. 
Wlien  constants  are  combined,  the  formula  for  drift  velocity  becomes: 


where 


w 

d 


2  -0.34 

V„  =  0.00925  V  d 
S  w 

surface  drift  velocity,  ft/sec 
windspeetl,  mph 
mean  depth,  ft 


(61 


Mean  surface  drift  velocities  are  shown  in  Table  4.  Arithmetic  mea^^;  of  drift 
velocities  were  used  for  analysis  of  wind  t>ffects  instead  (M  vi'ctor  sums, 
because  wind  tends  to  act  I’n  DMCAs  mi’re  as  a  mixing  agent  tlian  as  a  transp''rt- 
ing  agent.  An  exceptiim  would  he  ,i  strong  wind  blowing  toward  the  weir  of  a 


small,  shallow  area,  causing,  an  earlv  release  of  dve  at  h  i  g.h  concent  rat  i  mis 


r  Dve  Tracer  Tests 


Vector  suiTi  of  wind  measurements.  Directions  from  whicli  the  wind  blows,  measured  counterclockwise 
rom  site  axis.  Inflow  point  is  roughly  at  180°  and  outflow  is  at  0°. 

Storms  and  frequent  gusting  winds. 
l.VS  =  National  Weather  Service, 
p.c.  =  personal  communication. 


(Test  3  and  12).  As  Kif^ures  12  thrmigh  23  sliow,  m.uiv  of  tin*  tests  exper  ;  cu  c 
diurnal  cycling  of  wind  direction  tvpic.-il  <'1  coastal  areas.  Resultant  wind 

velocity  for  these  tests  is  quite  small,  but  wind  mixing  is  important. 

61.  Wind  conditions  varied  widelv,  with  average  velocities  ran>.’inc  Iro 

0  to  14  mph.  Most  resultant  wind  directions  tended  to  he  between  <)ii  and 

180  deg  relative  to  the  site  axes.  Temporal  variations  of  wind  speed'^  and 

directions  are  depicted  graphically  in  Figures  12  thri'ugh  23. 

Dve  Tests 

62.  In  each  test,  a  20-percent  solution  Hiv  weightl  Khodamine  W'l  f  luo- 
rescent  dye  in  water  was  used  as  a  tracer.  In  all  cases,  the  dve  was  added  t 
the  inflow  essentially  instantaneous  1 v  at  the  point  where  the  discharge  from 
the  dredge  pipeline  entered  the  basin.  The  amounts  used  ranged  between  I  and 
IS  gal,  but  was  usually  about  10  gal. 

63.  The  dye  concentration  at  the  overflow  weir  was  sometimes  measured 
hv  Dumping  water  continuously  through  either  a  Turner  Model  III  or  seriio 
10-005  fluorometer,  and  sometimes  hv  t.aking  grab  samples  at  irregular  inter¬ 
vals  and  measuring  the  dve  concentration  in  a  fluorometer  at  the  site  or  in 
the  lab  later.  In  one  case  fi’est  12,  black  Rock  Harbor),  both  lab  and  field 
measurements  were  made.  Results  were  reported  in  units  of  either  relative  civ 
concentration  (fluorometer  reading  or  cleflection)  or  absolute  dve  concentra¬ 
tion  (ppb).  The  raw  data  are  presented  in  digital  form  in  Appendix  ('. 

64.  The  duration  of  the  tests  ranged  1  rc'm  4  to  3'''f  hr,  but  most  were 
between  30  and  100  hr.  Test  durations  were  usually  between  one  and  two  times 
the  estimated  theoretical  volumetric  detention  time.  In  general,  test  dura¬ 
tions  were  not  sufficiently  long  enough  to  observe  decay  of  f  luoresience  to 
background  or  near-background  levels. 

Kffects  of  samp  1 e  hand  ling  and  dye  t  ype 

65.  The  effects  of  different  s.ample  handling  procedures  and  the  use  of 
different  dves  were  tested  by  running  replicate  experiments  at  the  I',  l.ick  Rock 
Harbor  site  fTest  I.').  Twc'  fluorescent  dves,  Rhc'damine  WT  and  1  issamine  i!  , 
were  s  imu  1  t  aneous  1  v  added  at  the  icilet.  S.impli-s  collected  at  t  lu'  outlet  weir 
W'ere  split  into  three  portions;  the  first  was  immediatelv  tested  tor  Kliod.i- 
mine  fluorescence,  and  the  other  two  were  returned  t  <>  the  I  a  bo  r  a  t  i  >  r .  Iliese 
second  and  third  peertions  of  each  sample  were  testec!  in  t  lie  laf-oialoi'.  sev'era 


weeks  after  collection  for  Rhodamine  and  lissamlne  fluorescence,  respectively. 
The  same  Turner  Series  lO-OOS  fluorometer  was  used  in  the  laboratory  as  in  the 
field.  Fluorescence  levels  were  read  in  the  laboratory  before  and  after  vig¬ 
orous  shaking  (2  to  3  sec)  to  determine  the  effect  of  resuspension  of  fine 
particles  on  fluorescence. 


PART  IV:  DATA  ANALYSIS  AND  KFSITTS 


("•verview 


66.  This  part  contains  a  description  of  the  procedures  that  were  u'ed 
to  process  and  reduce  the  residence  time  distribution  field  data  ccillected  as 
described  in  Part  III.  Results  of  the  experiment  conducted  at  the  Black  Pork 
site  to  determine  the  effects  <  f  dve  tvpe  and  sample  h.andlinp  are  repi^rted 
first.  Ihe  remainder  of  the  data  analvsis  had  as  its  objective  derivatic'n.  ('f 
a  method  for  predicting  hvdraulic  efficiencv  based  on  design  variables  sue'.  .)r 
site  geometrv,  discharge,  and  expected  or  "worst-case"  wind  conditions. 
Brieflv,  the  steps  followed  were: 

a.  Mf'st  of  the  dve  curves  were  incomplete  because  m.easur  em.en  t  ('1 
dve  co'ncent  rat  ion  was  not  continued  li'iig  enough  for  return  to 
background  levels  to  fe  observed.  The  missing  decav  linibs  or 
"tails"  of  the  dye  curves  were  estimated  b>.’  fitting  the  mea¬ 
sured  dve  concentration  values  to  a  log-normal  distribution 
curve  using  1  east -squares  best-fit  regression.  The  best-fit 
functions  were  then  used  to  predict  the  decav  limbs. 

h.  Next,  the  extrapolated  curves  were  used  to  calculate  numerical 
descriptors  such  as  mean  residence  times,  hydraulic  efficien¬ 
cies,  and  dispersion  indices. 

c.  A  matrix  of  correlation  coefficients  was  cu'mputed  for  the 
ntimerical  curve  descriptors. 

d.  The  mean  residence  time,  hvdraulic  efficiencv,  and  tiispersion 
index  were  used  as  dependent  variables  in  regression  analyses 
with  site  geometrv,  flow  rate,  and  wind  conditic'ns  as  indepen- 
rient  variables.  foefficients  were  computed  for  regression 
functions  based  on  both  dimensional  variables  and  dimensionless 
grc'ups  of  v.iriables. 

e.  Tiu'  '.’allies  of  mean  residence  time,  hvdraulic  efficiencv,  and 
dispersion  index  computed  using  the  regression  equatiiuis  were 
cor'.p.ired  with  observed  values.  ilvdraulie  efficiencies  and  dis¬ 
persion  indice'^  for  iivpotiiet  i  t  a  1  but  reasc'nabl  e  site  conditions 
were  (  "mp'ited  lining  tlie  regression  equations  and  were 

e  V , 1 1  M.i  t  ed  . 

A  simp'ifieii  predictive  equal  it>n  for  ti  vd  r  ,iu  1  i  (  effieienev  ,iv  ,i 
funiti'O'  ,if  iinlv  !  'W  w.is  fit  t  ('  .a  d.it.i  set  tinit  inoludi'd  tlie 

hita  .is  well  as  pui'Iistieii  d.it.i  for  ci-.I'ri'ie  (suit, ait  cfi.am- 
''•r  -,  w.isftw.itei  l.igoi'ns,  .ind  i  ph','‘  it. a!  mi'del. 
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Kt'focts  ot  Sample  Handling  and  Hve  Type 

61’.  Results  of  the  experiments  of  elfeots  of  dve  t'.-pe  and  sample  han¬ 
dling  are  sht’wn  in  Table  and  Fiptire  I’-i.  The  I.issam.ine  curves  produced  mt'an 
residence  times  and  dispersic'n  indices  that  were  rough!'.  'T)  percent  larger 
tlian  for  Rhodamine.  Tlie  decav  limbs  of  the  I.issamine  curves  liad  more  gradual 
slopes,  si.'gpestinp  tliat  the  l.is^aamine  did  not  degrade  under  the  ii’fluence  t'f 
ccmditions  in  tlie  DMtlA  as  cjuicklv  as  Rliodamine.  Slinking  the  samples  had 
insignificant,  though  I'bse rvah 1 e ,  effects  on  t  and  the  dispersion  index. 

11  i  t  f  erences  between  the  f  ;  e  1  il  and  lab  curves  were  also  slight  and  appear  to  be 
mainlv  because  of  a  higher  peak  concentration  measured  in  the  lab. 

(T  AK  VALUl-S 
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Figure  Residence  time  distributions 

mea' tired  simultaneously  using  two  dif¬ 
ferent  tracer  dves 

Fx  t  rapo  I  at  i  on  of  IncompK'te  Curves 

6R.  The  Id  dve  tracer  turves  were  classified  as  either  ci'-plete  .'r 
incomplete  (  i  ab  1  e  Ci.  tornlete  curves  were  those  that  !;ad  liec  o  ’  i''’"-  t':  it 
terminated  in  a  mild  !  I  uc  t  u.i  t  i  on  ahoiit  the  background  c.  I’i  c!  t  r,  t  i  '  '  .  'hie'-- 
groum!  concentrat  ioro,  determined  haseii  on  the  in-t  ia’  re  i  d  ;  e  i-.'-f 
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hr 

Background 

_ _ 

Max imum 

Concentration 

ppb 

So’ 

Di spersion 
I  ndex 

Rhodamlne  WT 

Field* 

8.0 

0.05 

53.45 

5.5 

0.82 

l.ab-unshaken 

7.0 

0.10 

199.90 

4.5 

0.99 

shaken 

6.0 

0.10 

179.90 

4.4 

1.01 

l.i  s?;am  ine 

l.ab-unshaken 

10.6 

20.0 

96.00 

9 . 2 

f).53 

shaken 

1  I  .0 

19.0 

83.00 

10.3 

0.52 

*  This  curve  was 

used  in 

the  analysis 

described  in  Tabl 

e  6  because 

all  other 

tests  were  run  with  Riiodamine. 

Table  6 

Dve  Curve  Analysis 


Coef  f ic i ent  s 


Complete  Regression  Coefficients  o f 

Test  or  Background  l  ^  Determination 

No.  Incomplete  Cone.  1  2  3  r- 


prior  to  dye  addition,  but  such  measurements,  if  any  were  made,  were  not 
reported.  Furthermore,  dye  concentrations  were  variously  reported  either  us 
relative  concentrations  (fluorometer  scale  readings)  or  in  absolute  concentra¬ 
tions  (ppb  of  dye),  thus  complicating  selection  of  reasonable  values  for 


background . 

69.  Most  curves  were  incomplete.  However,  since  they  exhibited  overall 
shapes  similar  to  the  log-normal  distribution,  log-normal  curves  were  fitted 
to  the  measured  points  using  a  regression  approach.  The  log-normal  distribu¬ 
tion  function  is  given  by 


C  =  - ! -  exp  < — ^  [ln(t)  - 

VTrr  ot  l2o  I 


(7) 


where  C  is  dye  concentration  and  t  is  time, 
constants,  this  expression  reduces  to 


If 


and 


are  taken  as 


ln(C)  =  b^  [ln(t)]^  +  ln(t)  + 


f«) 


which  is  a  quadratic  equation  for  ln(C)  in  terms  of  ln(t).  Values  for  the 
coefficients  b ,  ,  b 


2 


and  b  were  obtained  using  a  multiple  linear 

2 

regression  program  with  Inft)  and  [ln(t)]  as  the  independent  variables  and 

In(C)  as  the  dependent  variable.  Resultant  values  of  the  coefficients  b  , 

')  ^ 

b,  ,  and  b^  and  associated  coefficients  of  determination  r^  are  shown  in 
Table  6.  (’oef f  ic ients  were  determined  for  complete  as  well  as  incomplete 
curves  for  purposes  of  comparison.  Fitted  curves  and  measured  points  are 
shown  in  Figure  25.  The  fitted  curves  tended  to  give  rather  poor  fits  to  the 
measured  peak  dye  concentrations,  but  extremely  good  fits  to  the  decay  limbs; 
thus,  the  primary  objective  was  met.  The  regression  coefficients  were  used  to 
estimate  missing  values  for  decay  limbs  of  incomplete  curves.  Curves  were 
extended  either  to  background  concentration  or  three  times  the  theoretical 
residence  time,  whichever  occurred  first.  The  dye  curves  shown  in  Figures  12 
through  2i  include  estimated  values  for  Incomplete  curves.  Values  of  C 

o 

used  to  normalize  dye  concentrations  for  Figures  12  through  23  were  computed 
using  the  estimated  missing  values. 


TIME,  HR 


TIME.  HR 


TEST  3 


100 


TEST  4 


Numerical  Description  of  Dye  Curves 


Computation  of  descriptors 

70.  The  dye  tracer  curves  provide  representative  residence  time  distri¬ 
butions  from  the  tested  DMCAs .  Tn  order  to  quantitatively  describe  these 
distributions,  the  measured  times  and  concentrations  and  the  estimated  values 
for  missing  portions  of  incomplete  curves  were  entered  as  input  to  DYECON,  a 
dye  concentration  analysis  routine  of  the  Automated  Dredging  and  Disposal 
Alternatives  Management  System  (ADDAMS)  computer  program  CHayes  et  al.  1985). 
DYECON  computed  the  following  numerical  descriptors  for  each  curve: 


Mean  time 


mean  residence  time  ('Equation  2) 


Hydraulic  efficiency 


mean  residence  time  divided  by  the  theoretical 
residence  time  x  100% 


Minimum  time  =  the  time  at  which  the  first  dye  reading  above 
background  is  recorded 

Modal  time  =  the  time  at  which  the  maximum  observed  dye  reading  is 
recorded 


‘^lO’’^50’So 

Morrill  index 
Dispersion  index 


the  time  at  which  the  10,  50,  and  90  percentiles  of  the 
area  beneath  the  dye  curve  occur,  respectively 

So'^'^10 

dispersion  variance  as  computed  by  Equation  4 


DYECON  also  provides  a  variety  of  plot  options.  Full  documentation  for  the 
DYECON  program  is  found  in  Hayes  et  al .  (1985). 

71.  Results  of  the  DYECON  runs  are  in  Table  7.  Computed  mean  residence 
times  ranged  from  2  to  100  hr,  hydraulic  efficiencies  from  19  to  105  percent, 
and  dispersion  indices  from  0.03  to  0.60.  Hydraulic  efficiency  was  not  com¬ 
puted  for  Tests  9  and  II  because  pond  volume  (and  thus  T  )  was  unknown.  The 
two  impossibly  high  values  for  hydraulic  efficiency  (Tests  7  and  8,  Table  7) 
were  both  associated  with  curves  with  extremely  jagged  profiles  and  no  defi¬ 
nite  peak.  Although  practically  impossible,  those  two  values  were  retained 
and  used  in  subsequent  analysis  for  the  sake  of  consistency.  Marske  and  Boyle 
(1973)  also  reported  a  value  of  hydraulic  efficiency  in  excess  of  100  percent. 

72.  The  extremely  low  value  of  efficiency  for  Test  2  is  also  suspect. 
Estimates  of  site  volume  and  average  flow  rate  are  tenuous  for  this  test,  and 


Table  7 

Numerical  Curve  Descriptors 


Hydraulic 

Effi¬ 

ciency 

% 

Mini¬ 

mum 

Time 

hr 

Modal 

Time 

hr 

*^90 

hr 

Median 

Time 

hr 

^0 

hr 

Morrill 

1  ndex 

Di sper- 
s  ion 

I  ndex 

72.8 

12.2 

81.0 

28.4 

9.72 

8.3 

0.60 

ra 


the  normalized  dye  curve  has  a  shape  that  is  somewhat  atypical  relative  to  the 
others  (Figures  12  to  23). 

Comparison  of  descriptors 

73.  Relationships  among  dimensionless  forms  of  the  various  curve 

descriptors  (i.e.,  t/T  ,  t . /T  ,  t  /T  ,  etc.)  were  investigated  by  computing 

1  m 

a  matrix  of  linear  correlation  coefficients  (Table  8).  All  of  the  curve 

descriptors  except  t  .  /T  ,  the  Morrill  index,  and  the  dispersion  index  were 

min 

strongly  and  positively  correlated  with  t/T  .  The  Morrill  index  was  strongly 
and  positively  correlated  with  the  dispersion  index.  Additional  analysis  was 
therefore  focused  on  t  ,  t/T  ,  and  the  dispersion  index.  Marske  and  Boyle 
(1973)  recommended  use  of  the  dispersion  index  because  it  has  the  highest 
potential  for  statistical  reliability. 


Table  8 


Correlation 


Coefficients 


r  for  Numerical  Curve  Descriptors 


t 

t  . 
min 

'"modal 

*^50 

*^90 

'^10 

Mor¬ 

rill 

Disper- 
s  ion 

T 

T 

T 

T 

T 

T 

I  ndex 

I  ndex 

t/T 

1.00 

t  .  /T 
min 

0.21 

1  .00 

t  /T 

modal' 

0.58 

0.72 

1.00 

t  /T 

50' 

0.99 

0.32 

0.61 

1  .00 

t  /T 

90' 

0.97 

1 

o 

o 

o 

0.39 

0.93 

1  .00 

t  /T 
lO' 

0.63 

0.77 

0.98 

0.68 

0.43 

1  .00 

Morrill  index 

00 

o 

-0.54 

-0.69 

-0.32 

-0.12 

-0.68 

I  .00 

Dispersion  index 

-0.06 

-0.81 

-0.66 

-0.18 

0.15 

1 

o 

o 

0,84 

I  .00 

Dimensional  Variables 


74.  In  order  to  generate  a  predictive  relationship  for  use  in  DMCA 
design,  an  effort  was  made  to  study  the  relationship  of  mean  residence  time  to 
site  conditions.  The  mean  residence  time  of  water  in  a  DMCA  is  a  function  of 
geometry,  surface  drift  velocity  (produced  by  wind  shear),  the  wind  direction, 
and  the  mean  advective  flow  velocity  Q/WD  .  Algebraically, 

t  =  f(L,W,D,V  ,  ,  V  ,  cos  0)  (9) 

adv  s 

where 

t  =  mean  residence  time,  hr 

L,W,D  =  average  site  dimensions,  ft 

'^adv  ^  "’®3n  advective  flow  velocity,  Q/WD,  ft/sec 

Q  =  average  flow  rate  through  the  DMCA  during  the  dye  test  of 
period  of  Interest 

=  surface  drift  velocity  (liquation  b)  ,  ft/sec 

cos  f'  =  cosine  of  the  angle  between  the  site  axis  and  the  resultant 
wind  direction 

Tor  wind  blowing  directly  from  the  inlet  to  the  outlet  =  D  ,  and  for  wind 
blowing  directly  from  outlet  to  inlet  '  ■  =  IHO  .  u  is  measured  in  a 


counterclockwise  direction  from  the  site  axis  (a  line  from  inflow  to  outlet) 
to  the  direction  from  which  the  wind  blows. 


75.  A  correlation  matrix  was  computed  for  log-transf ormed  values  of  the 

variables  t  ,  d  ,  L  ,  W  ,  D  ,  V  ,  ,  V'  ,  and  cos  6  from  the  10  dye 

adv  s 

tracer  tests  for  which  theoretical  residence  times  were  available  (Table  9). 
The  value  of  zero  for  for  Test  8  was  transformed  by  assigning  a  value  of 

0.01  to  it  before  transformation.  Cos  6  was  transformed  by  using  the 
function: 


Transform  cos  6  =  sgn(cos  0)  *  ^log  [abs(cos  0)  +  .01 (10) 

because  6  of  90°  and  270°  should  have  the  same  effect  on  residence  time 
distribution. 


Table  9 

Correlation  Coefficients  r  for  Log-Transformed  Dimensional  Variables 


bog(t) 

Log(L) 

Log(W) 

Log(D) 

Log(V^) 

Log(cos  6)  ^'^^^'^adv^ 

Log(t) 

1.00 

Log(L) 

-0.002 

1 .00 

Log(W) 

-0.  10 

0.94 

I  .00 

Log(D) 

0.60 

0.59 

0.50 

1.00 

Log(V^) 

-0.08 

-0.15 

-0.20 

o 

1  .00 

Log(cos  0) 

0. 10 

0.04 

-0.10 

-0.05 

-0.19 

1  .00 

cn 

O 

-0.20 

O 

fs) 

-0.55 

0.16 

0.09  1.00 

76.  The  set  of  independent  variables  was  then  edited  to  eliminate 
redundant  or  collinear  combinations  (  r  greater  than  or  equal  to  0.5).  Tlie 
resultant  set  of  transformed,  uncorrelated,  independent  variables  (  W  , 

V^dv  ’  ''  ^  input  to  a  stepwise  regression  procedure  (Hunt  and 

Delagran  1984),  and  regression  coefficients  were  calculated  for  mean  residence 
time  as  a  function  of  the  independent  variables.  Ihe  resultant  equation  is: 

t  =  0.00019  r‘  -  o.s;  fin 

ndv 


A  correlation  matrix  for  log  transforms  of  these  variables  (Table  10)  revealed 
little  collinearity .  Nonpositive  values  were  log-transformed  as  before. 
Stepwise  regression  was  again  used  to  calculate  regression  coefficients  with 
the  results: 


r/T  .  0.29(v^/V^J-'’-"(l/u)»- 


94  2 

r  =  0.53 


and  dispersion  index: 


d  .  7.1 


r  =  0.26 


Wlien  the  Test  2  data  were  omitted,  stepwise  regression  yielded: 


/  .  No. 81/  .  No. 25/  N  -0.093  2  „ 

t/T  =  0.087^L/Wj  ('^s'^'^advj  ^ 


d  =  7.8 


(l/d)-" 


-0.53 


r  =0.24 


Table  10 

Correlation  Coefficients  r  for  Log-Transformed  Dimensionless  Variables 


LogCdispersion 

index) 

Log(t/T) 

Log(L/W) 

Log(L/D) 

Log(V  /V  ) 
s  adv 

LogCcos  0) 


Dispersion 

Index 


V  /V 

L/D  s'  adv 


cos  6 


78.  Predicted  and  observed  values  of  t/T  and  the  dispersion  index  are 
plotted  in  Figures  27  and  28.  Table  11  shows  observed  and  predicted  values 
for  Equations  12,  16,  and  17.  Figure  27  shows  that  Equation  16  predicted  t/T 
reasonably  well  for  all  tests  except  2  and  8.  Once  again,  the  extrapolation 


DISPERSION  INDEX 


Tab  1 e  II 

Observed  and  Predicted  Values 


Test 

Obs  t 

hr 

1‘red 

hr 

t* 

Obs  t/T 

I’red  t/T** 

Obs 

d 

Fred 

d 

1 

38.0 

1  7  . 

0 

7  3 

74 

0.60 

0.44 

O 

27.0 

49. 

0 

19 

66 

0.30 

0.27 

3 

1  .  7 

1  . 

3 

49 

62 

0.51 

0.52 

4 

5.9 

6. 

8 

63 

70 

0.27 

0.28 

5 

76.0 

54  . 

0 

49 

64 

0.60 

0.42 

6 

100.0 

163. 

0 

44 

53 

0.55 

0.30 

7 

43.0 

22  . 

0 

103 

103 

0.42 

0.49 

8 

1  72.0 

75. 

0 

105 

7  3 

0.55 

0.34 

9 

9.8 

— 

— 

— 

0.55 

— 

1  1 

41  .0 

— 

— 

— 

0.14 

— 

12 

8.0 

1  1  . 

0 

22 

21 

0.82 

0.83 

13 

7.0 

8. 

5 

51 

49 

0.08 

0.34 

*  Equation  12. 
**  Equation  16. 
Equation  17. 


Discussion  of  Results 

Reproducibility  of  dye  tests 

79.  The  observed  data  exhibit  considerable  scatter  about  values 
predicted  by  Equations  12,  16,  and  17  (Figures  26  through  2H  and  Table  1!). 
Mangelson  and  Watters  (1972)  observed  experimental  absolute  errors  ol  about 
10  percent  in  hydraulic  efficiency  while  working  with  an  indoor  model  basin 
40  ft  long  by  20  ft  wide  by  3.5  ft  deep.  There  was  no  wind,  and  sources  of 
error  were  more  tightlv  contrcilled  tlian  for  anv  of  the  ITUA  field  tests. 
Conversely,  Marske  and  B(')vle  (1973)  reported  that  the  results  nt  dve  tests  ii 
chlorine  ci.intact  ch.imbers  were  not  sensitive  t('  experimental  procedure,  am; 
"When  the  dve  test  was  repe.ated  on  the  same  basin  under  identical  1  low  condi¬ 
tions,  the  results  were  essentially  identic.il.”  Einht  of  the  jireilicted  vahu 
of  t/T  are  within  0 .  Ii  ol  the  observed  v.ilue,  with.  tw('  preilictions  (Jests  . 


and  8)  exhibiting  much  larger  errors  (0.46  and  0.34,  respectively).  Sources 
of  error  for  the  DMCA  tests  include  incomplete  data  for  decay  limbs,  the 
aforementioned  variation  in  methods  lor  determining  mean  depth,  inaccurate 
estimate  of  mean  flow  rates,  unsteadv  flows,  and  errors  introduced  by  using 
wind  data  collected  some  distance  awav  from  the  test  site. 

Predictions  for  extreme  conditions 

80.  Although  Fquations  13  and  16  provide  reasonably  good  fit  to  most  of 
the  observed  values,  thev  do  not  provide  realistic  values  when  values  of  inde¬ 
pendent  variables  outside  the  observed  range  are  used.  Although  Equation  12 
does  indicate  that  mean  residence  time  becomes  infinite  as  Q  approaches 
zero,  it  also  implies  that  mean  residence  time  is  completely  independent  of 
basin  length,  which  is  intuitivelv  false.  Equation  16  gives  impossibly  large 
hydraulic  efficiencies  when  even  fairly  modest  values  of  I./W  or  small  values 
of  D/W  are  input.  Hydraulic  efficiency  becomes  infinite  when  wind  speed  and 

thus  V  are  zero, 
s 

81.  Since  Equation  16  is  based  on  dimensionless  variables,  predictions 
from  it  can  be  compared  with  data  from  other  studies.  Unfortunately,  pub¬ 
lished  results  were  found  for  only  two  other  similar  investigations;  a  study 
of  chlorine  contact  chambers  by  Marske  and  Boyle  (1973)  and  a  study  of  model 
and  prototype  waste  stabilization  ponds  by  Mangelson  and  Watters  (1972). 

Marske  and  Boyle  include  wind  and  depth  observations  for  only  two  of  their 
tests.  Most  of  Mangelson  and  Watters's  data  were  from  an  indoor  model  pond 
with  no  wind.  Although  they  did  no^e  the  presence  of  wind  effects  for  their 
prototype  tests,  they  did  not  record  wind  speeds. 

82.  Marske  and  Boyle  (1973)  reported  values  of  hydraulic  efficiency  of 

83  and  84  percent  for  tests  run  in  a  basin  336  ft  long  by  16  ft  wide  by  3  ft 

deep.  Details  of  the  two  tests  were  virtually  identical  except  for  wind 

direction.  Wind  was  20  mph  directly  downstream  in  the  first  test  and  20  mph 

directly  upstream  for  the  second.  Equation  16  predicts  a  hydraulic  efficiency 

of  280  percent.  The  ratios  I./W  and  I)/W  were  both  an  order  of  magnitude 

greater  for  these  tests  than  for  the  tests  on  DMOAs,  and  the  ratio  \’  /V  , 

s  a(i\' 

was  smaller  than  for  all  DMCA  tests  except  the  one  with  zero  wind  speed 
(Test  7).  It  theretore  appears  that  E(|uation  13  is  not  generally  applicabU’ 
outside  the  range  ot  conditions  enci)untered  in  the  DMt'A  tests,  even  though  it 
is  based  on  dimensionless  variables. 


83.  Equation  12  is  heavily  intluenced  by  width  and  is  insensitive  tc' 
changes  in  wind  speed.  I’se  of  Equation  12  to  estimate  tlie  et  tents  nt  adding 
spur  dikes  is  particularly  misleading  because  length-width  ratiii  is  not 
included  in  any  form.  Similarly,  Equation  lb  gives  unreasonable  values  for 
hydraulic  efficiency  when  length-width  ratios  exceed  maxirami  observed  values 
(3  to  4) . 

84.  In  summary,  the  regression  Equations  12  and  lb  give  unreasonable 

values  where  used  outside  the  observed  ranges  of  the  independent  variables. 

However,  they  do  highlight  the  influence  of  L/W  on  DMCA  mean  residence  time 

and  hydraulic  efficiency.  (In  Equation  12,  note  that  V  ,  =  (i/U’D  .)  They 

adv 

also  highlight  the  lack  of  influence  exerted  bv  other  variables  such  as  wind 
speed  and  direction  on  hydraulic  efficiency.  Less  definitive  results  were 
obtained  for  the  dispersion  index. 

Relative  Effects  of  Site  Condition  Variables 

Effect  of  length-width  ratio 

85.  Regression  equations  based  on  only  DMCA  data  fail  to  adequately 
describe  effects  of  changing  L/W  on  hydraulic  efficiency  because  the 
observed  range  of  L/W  was  so  small.  The  m.odel  study  by  Mangelsnn  and 
Watters  (1972)  and  the  chlorine  contact  chamber  study  by  Marske  and  Boyle 
(1973)  made  measurements  of  t/T  over  a  much  wider  range  of  L/W  .  Both  sets 
of  investigators  found  a  direct  relationship  between  L/W  and  hydraulic 
efficiency.  Marske  and  Boyle  (  1973)  recommended  use  of  I./V,’  =  40  for  design 
of  chlorine  contact  chambers.  Using  regression  equations  based  on  tlieir  data, 
they  predicted  a  hydraulic  efficiency  of  81  +  24  percent  and  a  dispersion 
index  of  0.02  for  such  a  basin.* 

8b.  Figure  29  shows  all  of  the  available  data  for  the  effect  of  L/W 
on  hydraulic  efficiency.  Data  from  UMCAs,  chlorine  contact  chambers,  waste- 
water  lagoons,  and  the  physical  model  described  previously  are  included. 
Regression  coefficients  were  calculated  using  these  data  for  the  following 
equat ion : 


*  The  dispersion  index  wyis  estimated  from  a  plot  of  d  vs.  I./W'  .  Hydraulic 
efficiency  was  predicted  using  a  regression  c<]uation  with  d  as  tlie 
independent  variable  and  t/T  as  the  dependent  variable. 


is  proportional  to  the  -0.3  and  -0.25  power  of  depth,  respectively.  This  con¬ 
firms  the  findings  of  Mangelson  and  Watters  (1972).  In  other  words,  increases 
in  depth  for  a  given  width  tend  to  increase  the  fraction  of  volume  in  dead 
zones . 

88.  On  the  other  hand,  depths  that  are  too  shallow  can  degrade  DMCA 
performance.  Montgomery  (1978)  suggests  2  ft  as  a  minimum  DMCA  ponding  depth 
and  attributes  some  observed  DMCA  short-circuiting  to  inadequate  ponding 
depths.  Walski  and  Schroeder  (1978)  present  a  design  procedure  for  DMCA  out¬ 
let  weirs  that  links  effluent  quality,  weir  length,  and  ponding  depth.  There¬ 
fore,  DMCA  depths  should  generally  be  2  to  8  ft. 

Effect  of  wind  speed  and  direction 

89.  In  Gallagher  and  Company  (1978),  the  results  of  two  DMCA  dye  tests 

in  the  same  DMCA  indicated  that  a  10-  tn  I  S-mnb  wind  tov.’ard  the  v.’eir  reduced 

t^/T  by  0.13  relative  to  a  no-wind  condition.  However,  the  dye  "curve"  for 

the  wind  test  was  only  two  points.  Marske  and  Boyle  (1973)  concluded  that 

wind  effects  on  performance  of  a  basin  with  L/W  =  21  were  relatively 

unimportant  based  on  the  results  described  above.  Although  the  20-mph  wind 

toward  the  weir  resulted  in  a  reduction  of  dimensionless  initial  time  t./T 

1 

of  0.18  relative  to  a  test  with  20-mph  wind  away  from  the  weir,  the  other  dye 
curve  descriptors  were  unchanged.  Conversely,  Mangelson  and  Watters  (1972) 
noted  that  simultaneous  dye  tests  in  two  adjacent  waste  stabilization  ponds, 
with  one  oriented  transverse  and  the  other  parallel  to  the  wind  produced  sig¬ 
nificantly  different  results.  The  pond  parallel  to  the  wind  (wind  blowing 
from  inlet  to  outlet)  experienced  "serious  short-circuiting."  Quantitative 
results  were  not  reported.  Data  for  the  DMCAs  were  insensitive  to  wind  direc¬ 
tion,  although  the  high,  early  spikes  in  curves  for  lests  3  and  12  are  doubt¬ 
less  the  result  of  wind  blowing  across  relativelv  small  areas  toward  the  weir. 
[•:<luation  16  indicates  tliat  hvdraulic  elticiencv  is  proportional  to  wind  speed 
to  the  -(1.2  power.  Based  on  the  above  evidence,  it  seems  tliat  wind  effects  on 
hydraulic  elliciencv  are  slight,  but  wind  effects  on  the  overall  residence 
time  distribution  are  presently  uncertain,  but  probable  s i  >'n i f  i  cant  . 

Effect  of  inlet  and  outlet  configuration 

90.  Mangelson  and  Watters  ('!<*;■;’')  found  model  pi'iul  pi' r  i  o  ritianc  e  most 
efficient  when  '.'ertical  or  horizontal  diffusers  placed  at  periini'ter  locations 
were  used  to  distribute  inflow.  Most  chan>’,es  in  t  'T  resulting  froia  cli.anges 

u  t  i  i' t  c  c'l;  1  i  >;u  r.i  t  i  ons  w/ere  O.I  or  less,  althongb,  the 


in  inlet  .ind  /or 


conf igurat ion  most  similar  to  a  typical  rectangular  DMCA  had  the  lowest  value 
of  t/T  (0.456),  and  one  configuration  using  a  long  diffuser  for  inflow  had 
t/T  =  0.672.  Marske  and  Boyle  (1973)  demonstrated  th  at  a  sharp-crested  weir 
is  superior  to  a  Cipolleti  weir  and  recommended  minimization  of  weir  overflow 
rates  for  greater  hydraulic  efficiency.  The  DMCA  data  base  did  not  contain 
sufficient  data  to  adequately  compare  the  effects  of  various  types  of  inlets 
and  outlets.  Walski  and  Schroeder  (1978)  present  guidance  for  the  design  of 
DMCA  weirs  based  on  typical  solids  concentration  profiles  and  weir  withdrawal 


zones . 


PART  V:  DESIGN  GUIDANCE 


I ntroduction 


91.  Residence  time  distribution  is  an  important  consideration  in  the 
design  of  most  DMCAs.  All  but  the  largest  and  most  permanent  DMCAs  experience 
periods  when  residence  time  distribution  is  critical  to  effluent  quality,  usu¬ 
ally  near  the  end  of  filling,  when  DMCA  pond  volume  and  surface  area  are 
least.  If  a  DMCA  is  sized  large  enough  to  provide  adequate  effluent  quality 
even  when  hydraulic  efficiency  is  low  and  dispersion  is  high,  consideration  of 
residence  time  distribution  in  design  is  not  necessary.  However,  such  a  DMCA 
design  is  economically  inefficient  with  regard  to  effluent  quality  constraints 
on  size.  As  shown  in  Part  1,  improving  hydraulic  efficlencv  from  40  to 

60  percent  reduces  the  required  ponded  volume  by  33  percent. 

92.  Plug  flow  represents  the  ideal  residence  DMCA  time  distribution. 
Although  this  ideal  is  not  achievable  in  practice,  low  dispersion  indices  and 
high  hydraulic  efficiency  should  be  the  designer's  goal.  Design  decisions 
generally  determine  residence  time  distribution;  operational  actions  to 
improve  effluent  quality  such  as  reducing  flow  rates  or  increasing  ponding 
depth  are  often  infeasible  because  of  other  constraints. 

Determining  Residence  Time  Distributions  for  Existing  DMCAs 

93.  The  residence  time  distribution  of  an  existing  DMCA  may  be  deter¬ 
mined  by  performing  a  dye  tracer  test  during  disposal  operations  with  flow 
rates,  wind  conditions,  ponding  depths,  and  surface  areas  similar  to  the 
conditions  of  interest.  Recommended  procedures  for  dye  tracer  tests  are  pro¬ 
vided  in  Appendix  B  of  this  report.  If  the  results  of  such  tests  indicate 
inadequate  hydraulic  efficiency,  residence  times  may  be  increased  by  reducing 
inflow  rates  or  by  using  the  tecliniques  for  improving  hydraulic  efficiencv. 

94.  In  the  absence  of  site-specific  field  data,  the  liydraulic  effi¬ 
ciency  of  a  given  site  may  be  estimated  using  Equation  19  from  Part  IV. 


t/T  =  0.9  [l 


exp  (-0.3  L/W)J 


(191 


This  relationship  is  plotted  in  Figure  30  along  with  the  DMCA  data.  Values  of 
t/T  from  this  equation  should  be  rounded  to  the  nearest  O.I.  L  for  use  in 
this  equation  is  the  distance  from  the  DMCA  inlet  to  the  outlet  measured  along 
the  center  line  of  the  flow,  and  W  is  the  average  width  of  the  flow  path. 


Figure  30.  DMCA  design  curve.  Hydraulic  efficiency  as  a 
function  of  I  'W  .  Solid  curve  is  Equation  19 

95.  The  hydraulic  efficiency  correction  factor  HECF  will  be  the  recip¬ 
rocal  of  t/T  .  If  the  correction  factor  is  used  in  additional  analysis  for 
the  prediction  of  effluent  quality,  the  engineer  should  be  careful  to  multiply 
the  HECF  by  the  DMCA  pond  volume  (or  surface  area  for  zone  settling  analysis) 
likely  to  exist  at  the  time  of  interest  instead  of  its  initial  condition.  For 
example,  use  of  the  surface  area  enclosed  by  the  perimeter  dikes  could  intro¬ 
duce  error  since  DMCA  surface  area  is  usually  significantly  reduced  during 
filling  by  the  formation  of  deltas  or  fans  of  coarse  material  at  Llie  points  of 
inflow. 


Predicting  Hydraulic  Efficiency  of  Proposed  DMCAs 

96.  Hydraulic  efficiency  is  incorporated  into  sizing  o|  pi'oposed  DMiVXs 
by  multiplying  the  required  surface  area  fcir  effective  si'ttling  (as  deteniuned 
from  lab  tests)  by  the  reciprocal  of  livdraulic  efiicienc',,  il  -one  •  ^  t - 

tling  analysis  is  employed,  or  bv  multiplvie.-  t!ie  t.'sniii.  ii'  liiv.i  t  i . 
the  HECF  if  f  1  occu  1  en  t  se  1 1  1  i  ng  .  ina  1  ■  i  s  i  '  i  ■  .  c  ,  ;  '  i  s  a  ,  ,  .  .  .  i  i  ,  i 


mal  ly  exhibit  zone  settling,  and  that  f  ■;  i- 


where  saltwater  sediments  are  to  be  dredged.  A  detailed  description  is  given 
by  Palermo,  Montgomery,  and  Poindexter  (1978).  The  floccnlent  settling  analy¬ 
sis  is  normally  used  for  design  of  DMCAs  to  receive  freshwater  sediments  with 
significant  clay  fraction  and  for  situations  where  sediments  exhibit  zone 
settling  but  DMCA  effluent  contaminant  concentrations  are  of  concern.  A 
detailed  description  of  the  procedure  and  an  example  are  given  by  Palermo 
(1985)  . 

97.  Thus  in  the  zone  settling  analysis,  the  design  surface  area  A^  is 
equal  to: 


A  =  HF.CF  *  A,  =  (  — ]  *  A  (20) 

'■  Vtn/ 

where  A^  is  the  required  surface  area  determined  from  laboratory  tests.  For 

flocculent  settling,  the  design  residence  time  T^  is  equal  to  the  required 

d 

residence  time  determined  from  laboratory  tests  T^^  multiplied  by  the  HECF : 

T_,=  HECF  *  T  =  (—]  *  T„  (21) 

98.  A  simple  approach  for  evaluating  the  layout  of  a  proposed  DMCA 

would  be  simply  determining  L/W  from  site  plans  and  using  Equation  19  to 
predict  t/T  .  As  noted  above,  I,  should  be  the  distance  from  inlet  to 
outlet  along  the  center  line  of  the  flow,  and  W  should  be  the  average  width 
of  the  flow  path.  The  computed  value  of  t/T  could  then  be  used  in  either 
Equation  20  or  21  (whichever  is  appropriate)  to  determine  if  the  DMCA  size  is 
adequate  to  ensure  good  effluent  quality.  Design  would  proceed  by  trial  and 
error:  (a)  choose  a  layout,  (b)  compute  efficiency,  (c)  check  site  area  or 

volume,  and  (d)  repeat  if  necessary.  The  DMCA  area  (or  volume  and  flow  rate) 
used  in  Equation  20  or  21  should  be  values  representative  of  critical 
conditions. 

99.  However,  since  hydraulic  efficiency  and  surface  area  of  square  or 
rectangular  DMCAs  are  both  functions  of  I-  and  W  ,  a  more  appropriate  design 
procedure  would  be  to  combine  Equations  18  and  19  or  20  to  economically  opti¬ 
mize  DMCA  dimensions  and  the  number  of  spur  dikes.  A  simplified  procedure  for 
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such  a  cost  analysis  is  given  in  the  following  section,  A  similar  analysis 
might  be  used  to  evaluate  the  economics  of  adding  spur  dikes  or  baffles  to  an 
existing  DMCA. 


Improving  Residence  Time  Distributions 

100.  Site  conditions  that  promote  nonuniform  flow  tend  to  degrade  resi¬ 
dence  time  distribution  away  from  the  plug-flow  ideal.  To  promote  uniform 
flow,  the  bottom  of  the  DMCA  should  be  free  of  vegetation  larger  than  grasses. 
Major  topographic  variations  such  as  ridges  and  swales  should  be  graded 
smooth.  Inlets  and  outlets  should  be  designed  and  located  so  as  to  avoid  con¬ 
centration  of  flow  and  adjacent  dead  zones. 

Spur  dikes 

101.  Baffles  or  internal  spur  dikes  are  by  far  the  most  effective  and 
reliable  devices  that  can  be  used  to  significantly  reduce  short-circuiting  and 
dispersion.  One  well-placed  spur  dike  can  nearly  double  the  length  and  halve 
the  width,  thereby  increasing  the  length-to-width  ratio  by  a  factor  of  3  to  4. 
Figure  31  shows  several  good  spur  dike  arrangements.  In  general,  spur  dike 
length  should  be  0.75  times  the  length  of  the  side  parallel  to  the  dike.  Spur 
dikes  should  subdivide  the  DMCA  pond  into  sections  with  near-equal  surface 
areas.  Equation  19  indicates  that  L/W  ratios  above  10  provide  very  little 
improvement  in  the  hydraulic  efficiency,  and  therefore  use  of  spur  dikes  to 
achieve  L/W  ratios  greater  than  10  is  probably  unnecessary.  Extremely  high 
L/W  ratios  might  produce  advective  velocities  high  enough  to  scour  and  resus¬ 
pend  deposited  sediments. 

102.  Spur  dikes  force  the  flow  to  reverse  direction,  which  can  have 
several  advantages.  The  first  turn  breaks  the  momentum  caused  by  the  inflow 
jet  and/or  wind,  and  each  subsequent  turn  redistributes  the  flow,  eliminating 
any  tendency  for  the  flow  to  concentrate  in  a  narrow  hand.  The  influence  of 
wind  on  short-circuiting  is  reduced,  since  the  wind  cannot  establish  a  surface 
current  that  travels  directly  from  the  inlet  to  the  outlet.  The  wind  is 
forced  to  blow  across  or  against  the  flow  over  a  portion  ol  the  Mow  path, 
which  reduces  large  wind-induced  circulation  patterns. 

103.  Commercially  available  floating  b.iffles  may  oiler  several  advan¬ 
tages  to  permanent  spur  dikes  (ligure  321.  Ihev  can  be  installed  quiiklv, 
even  after  disposal  has  begun.  Thev  can  be  moved  after  deltas  tonii  to 


a.  Longitudinal  Spur  Dikes 


b.  Transverse  Spur  Dikes 


c.  Combination  Spur  Dikes  for  Difficult 
Access  and  Shape  Restrictions 


Figure  31.  Examples  of  longitudinal  and  transverse  spur  >iike  rontii’- 
urations.  M  is  F/W  without  spur  dikes  and  M*  approximate 
I,/W  ratio  with  spur  dikes  (from  Callagher  and  Comiianv  (Ih'H'l't 


increase  the  L/W  ratio  near  the  end  of  the  disposal  life  of  the  DMCA. 

Unlike  spur  dikes,  baffles  use  very  little  of  the  DMCA  volume,  thus  allowing 
more  volume  for  storage  and  settling. 

Inlet  devices 

104.  The  most  important  design  aspect  of  inlet  devices  is  the  necessitv 
of  destroying  the  forward  momentum  let  and  distributing  the  flow'  throughout 
the  width  of  the  basin.  Ways  to  achieve  this  design  goal  include  the 
foil  owing : 

a.  Attach  a  flow  splitter  to  the  end  of  the  dredge  pipe  to  split 
the  flow  and  direct  it  toward  the  sides  of  the  DMCA. 

b.  Attach  a  spoon  to  the  end  of  the  dredge  pipe  to  spread  the  flow 
as  a  fan  in  all  directions  from  the  inlet. 

£.  Use  a  wye  or  a  manifold  to  distribute  the  inflow  across  the 

entire  inlet  end;  preferably,  each  point  of  inflow  should  have 
a  device  to  break  the  forward  momentum. 

d.  Install  a  splash  baffle  in  front  of  the  inlet  pipe. 

Inlet  devices  like  manifolds  that  divide  the  flow  should  be  designed  to  avoid 
clogging.  During  many  disposal  operations,  it  is  necessary  to  move  the  dredge 
pipe  to  different  inflow  points.  In  such  cases,  special  inlet  devices  must  be 
portable . 

Outlet  devices 

105.  The  ouvjlet  device  for  a  DMCA  is  typically  a  weir.  Weirs  concen¬ 
trate  flow  In  the  area  directly  in  front  of  them,  creating  dead  ?,ones  on 
either  side  (Figure  33).  The  size  of  the  dead  zones  can  be  reduced  bv 
increasing  the  length  of  the  weir  or  by  installing  several  weirs  along  the 
outlet  end.  However,  care  must  be  taken  to  ensure  that  all  of  the  w'eir  crests 
are  at  the  same  elevation  to  avoid  concentrating  the  flow  at  locations  where 
the  weir  crests  are  lower.  Walski  and  Schroeder  (I97R)  provide  guidance  and 
design  nomograms  for  DMCA  weir  design. 

Placement  of  inlets, 
outlets,  and  spur  dikes 

106.  The  objective  of  the  placement  of  inlets,  outlets,  and  spur  dikes 
is  to  reduce  wind  effects,  maximize  the  I./W  r£itio,  aiul  ensure  that  water 
flows  through  the  entire  DMCA  pond,  keeping  dead  zones  to  a  minimum.  I'iie 
placement  of  the  inlet  and  outlet  is  frequentlv  controlled  bv  site  and  project 
constraint.s  such  as  the  location  with  respect  to  the  eiredge  anil  reci‘ivinr 
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water  body.  Nevertheless,  efforts  should  be  made  to  maximize  the  distance 
between  the  inlet  and  outlet  and  to  place  the  Inlet  and  outlet  on  opposite 
sides  of  the  DMCA,  unless  they  are  separated  by  a  spur  dike. 


Cost  Analysis 

107.  The  surface  area  and  length-width  ratio  for  a  square  or  rectangu¬ 

lar  DMCA  are  both  functions  of  length  and  width.  If  baffles  or  spur  dikes  are 
installed,  the  number  and  dimensions  of  the  spur  dikes  enter  the  functions, 
but  area  and  L/W  are  still  functionally  related.  By  means  of  equations  to 
describe  these  relationships,  the  least-cost  configuration  of  site  length, 
width,  and  number  of  spur  dikes  may  be  determined  for  a  given  value  of  site 
area  A^^  (or  for  f  locculent  analysis)  . 

Constraint  function 

108.  The  following  relationships  may  be  used  to  transform  Equations  20 
and  21  into  functional  relationships  between  1,  ,  W  ,  and  the  number  of  spur 
dikes  or  baffles  N  .  Assuming  a  rectangular  DMCA  with  points  of  inflow  and 
outflow  at  the  midpoints  of  the  shorter  sides  (Figure  34),  pond  surface  area  = 
pond  length  x  pond  width  -  area  occupied  by  spur  dikes,  or 

A,  =  LW  -  I.’NW  L  (22) 

d  s 


where 

L  =  pond  length,  inlet  to  outlet,  ft 
W  =  average  pond  width,  ft 

N  =  number  of  spur  dikes  (all  spur  di.  !s  assumed  to  be  longitudinal) 

W  =  width  of  spur  dikes  at  water  surface  elevation 
s 

The  theorectical  residence  time  is 

T  =  volume/Q  (231 

d 

where 

ft  =  average  flow  rate,  cfs 

Pond  volume  is  porn!  length  “  pond  width  '  pond  depth  -  volume  occupied  bv  spur 
dikes,  or 


Pond  volume  =  I.U'D  -  P'N\’  I 


(2  41 


VSVN 


HFCF  =0.9  1  -  exp  -0. 3L ' (L/W) (N  +  I)' 


(^9) 


By  substituting  Kquation  25  into  Equation  20,  the  following  is  derived: 


.(W  -  L'NW  ) 
s 


0.9 


1  -  exp 


-0.31.'  (L/W)  (\  +  1)' 


Substituting  Equation  25  into  Equation  21  gives: 


L(WI)  -  L'N'V  ) 
u 


0.9 


1  -  exp 


-0.3L' (L/W) (N  +  1)' 


L(W  -  aN)  = 


1  -  exp  c(L/W)(K  +  !)■ 


where  a  ,  b  ,  and  c  are  constants.  If  W>>aN,  then 


(26) 


(27) 


I  10.  If  all  variables  in  Equations  26  and  27  are  fixed  except  N'  ,  I.  , 
and  W  ,  the  two  equations  become  identical: 


(28) 


LW  " 


I  -  exp 


_b _ 

c(L/W)(N  +  I)“ 


(29) 


Equations  28  and  29  mav  be  referred  to  as  constraint  functions,  since  OMCA 
design  must  satisfy  tliese  equations  in  order  to  provide  acceptable  effluent 
qua  1 i ty . 

Lost  function 

Ill.  The  cost  of  a  OMf’A  consists  of  the  cost  of  dikes,  land,  and  other 
items  like  outlet  structures,  channels,  access  roads,  and  miscellaneous  items 
The  area  of  land  required  (excluding  land  necessarv  for  access  and  rights-of- 
wav)  will  be  (Figure  I'll: 


7  i  t  e  are  a 


(I.  +  :’w)(w  +  :'w) 


(  7(1'! 
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Figure  iS.  Definitiun  sketcTi  i Hr  cost  analvsis 

where  w  is  the  base  width  of  tfie  perimeter  dike. 

The  lengtii  of  the  perimeter  dike  will  he  (Figure  dSl: 

Perimeter  dike  length  =  J(I.  ■*  .'w  +  Wl 

The  total  length  of  the  spur  dikes  will  h»'  (Figurt'  j  : 


Spur  dike  length  -  NT  'I. 

If  1  is  the  ratio  of  tlie  unit  cost  of  the  spni  dike  to.  'in- 
perimeter  dike,  then  the  cost  of  the  hVt'A  .likes  and  i.irsl  will 


$ 


1^2  (L  +  2w  +  W)  +  fNL'LjC^ 


[(L  +  2w)(W  +  2w)J 


(33) 


where  C  is  the  cost  of  the  perimeter  dike  in  $/ft  and  C  is  the  cost  of 
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land  in  $/ft  . 

112.  A  designer  may  optimize  the  portion  of  DMCA  costs  represented  by 
Equation  33  subject  to  the  constraint  imposed  by  Equation  28.  In  other  words 
the  least  costly  combination  of  L  ,  W  ,  N  ,  and  L'  that  will  provide 
acceptable  effluent  quality  may  be  selected.  Some  simplified  examples  follow 
The  basic  assumptions  in  these  examples  are  the  same  as  those  used  by 
Gallagher  and  Company  (1978)  in  their  Appendix  D. 

Basic  assumptions 

113.  The  following  assumptions  are  used  in  all  examples: 

a.  All  spur  dikes  are  longitudinal,  and  L'  ,  the  ratio  of  spur 
dike  length  to  pond  length,  is  0.75. 

b.  Perimeter  dikes  have  side  slopes  of  1:3,  crown  widths  of 

~  3 

10  ft,  and  heights  of  7  ft.  Their  unit  volume  is  217  ft  /ft. 

c.  Spur  dikes  have  side  slopes  of  1:2,  crown  widths  of  10  ft,  an 
heights  of  7  ft.  Spur  dike  width  at  the  water  surface  eleva¬ 
tion  is  10  ft,  and  unit  volume  below  the  waterline  V  is 

3 

100  ft  /ft, 

d.  Ponding  depth  D  is  fixed  at  5  ft,  so  all  dikes  have  2  ft  of 
f  reeboard , 

e.  The  average  flow  rate  Q  is  27  cfs. 

f.  The  required  residence  time  from  laboratorv  tests  (floe- 

— 

culent  analysis)  is  44.75  hr,  while  the  required  surface  area 
A.  from  lab  tests  (zone  settling  analvsis)  is  20  acres. 

The  unit  cost  of  perimeter  dikes  C^  is  $Ib/ft.  Unit  cost 

for  spur  dikes  is  half  as  much,  so  f  =  0.5  . 

h.  The  unit  cost  of  land  is  $500/acre  or  ?0.011/ft  . 

i.  The  DMCA  is  rectangular  in  shape  with  uniform  depth  so  that 

A  =  TU'  and  volume  =  I.KD  . 


l-X.imp  I  e  l--neitlier  length  nor  width  fixed 

114.  In  the  case  that  neither  DMCA  width 
site  cundltioii'  ,  costs  .are  a  1  unction  nt  ! 

V, lines,  Iqnaf  I  ni  r,a','  lie  used  ti-  ci'mtuite  I'.-.' 

with  the  ri'.'en  \,alue  of  1  /U'  to  !  ind  !  and 
lo  ed  t  i'  11  "d 


nor  lencth  are  e on ,s t  r ,i  i  ned  '' 
and  N  .  (  iven  N  and  ' 

.  I,'.-.'  IT', a'.’  then  l-e  I'omhii'ed 

V.'  ,  .'ind.  then  '  i  ;\!.i  t  1  .mi  '  ^  ■  i: 

and  r  1  c n  r  e  ■  sh ow 


the  total  cost  of  1  and  and  d  1  lei  . 


1  .ih  1  t 


the  effects  of  L/W  and  N  on  the  total  cost  of  land  and  dikes  under  the  set 
of  assumptions  listed  in  paragraph  113.  Both  Table  12  and  Figure  36  show  that 
spur  dikes  reduce  costs  when  the  above  assumptions  hold.  However,  the  incre¬ 
mental  rate  of  return  on  spur  dikes  diminishes  rapidly,  so  one  to  three  longi¬ 
tudinal  spur  dikes  will  be  the  best  arrangement  for  most  situations. 

Table  12 

Minimum  Cost  Configuration  for  DMCA  with  Neither  Length 


«  100 
o 

(J 


N  IS  NO,  OF  SPUR 
DIKES. 


?.5  5.0  7.5 

POIMH  L  W  wnnoui  SI'UH  lUKtS 


Figure  36.  Relationship  of  cost  of  DMCA  land  and 
dikes  to  the  number  of  spur  dikes  and  tlie  overall 
nond  L/W  ratio 


Example  2 — width  fixed 

115.  In  most  cases,  the  designer  cannot  freely  vary  the  length  and  width 
of  the  DMCA  because  available  tracts  of  land  are  limited  in  size  and  shape. 

In  a  case  where  the  pond  width  of  the  proposed  DMCA  is  set  at  300  ft  but  the 
length  may  take  on  any  value  greater  than  or  equal  to  300  ft,  Equation  29  may 
be  solved  by  trial  and  error  for  L  using  several  different  values  of  K  , 
the  number  of  spur  dikes.  Table  13  shows  the  effect  of  N  on  DMCA  cost.  A 
comparison  of  Tables  12  and  13  shows  that  restricting  DMCA  shape  can  be 
costly . 


Table  13 

Variation  of  Cost  with  Number  of  Spur  Dikes  for  DMCA 
with  Width  Fixed  at  300  Ft 


Spur  Dikes 

Overall 

No. 

L/W 

Cost,  $K 

L,  ft 

W,  ft 

LW,  acres 

.7 

140.9 

3,500 

3 

.0 

153.4 

3,300 

3 

.3 

178.0 

3,400 

3 

.7 

203.9 

3,500 

3 

IM 
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116.  In  summary,  the  relationship  between  L/W  and  hydraulic 
efficiency  allows  economic  optimization  of  DMCA  layout — pond  shape,  size,  and 
the  number  of  spur  dikes.  Floating  baffles  are  often  economically  attractive 
relative  to  earthen  spur  dikes  both  because  of  their  cost  and  because  they 
ocr  ,py  so  little  of  the  DMCA  pond  volume.  Certainly  more  sophisticated 
approaches  to  optimization  with  fewer  simplifying  assumptions  (such  as  fixing 
pond  depth)  are  warranted  for  some  projects. 


PART  V  1  ;  CONf 


117.  Flow  residence  time  <1  i  st  r  :  ‘  ni  : 

DMCA  design.  Underestimation  of  t'.vdr.in!  ii  i  •  ‘  \ 
increase  cost,  whereas  overestimation  co-n  rc-v''  ■  ■ 

effluent  quality.  Current  design  pri'ccciuM-  •  o  '  ;  i  r 

ef  f  iciency  instead  of  the  overnl  1  residence  t  ime  .i  :  •  •  r  .  •  "  ■  .  ■  ,  ,  , 

quality  is  determined  by  the  product  of  tt;e  soli, Is  r, ' 
and  the  residence  time  distribution. 

118.  DMCAs  tend  to  have  log-normal  residence  time  d  i  <  r  r  : 'ot  : -r  ■  .  ■  f 

DMCAs  have  hydraulic  efficiencies  between  20  and  '()  percent  witf  ,!  i  spo  c  i  <  c-; 

indices  between  0.3  and  0.6.  Wind  tends  to  affect  the  minirium  and  rmsla!  resi¬ 
dence  times,  and  thus  the  overall  shape  of  the  residence  time  d i st r ilnit ion , 
but  not  the  hydraulic  efficiency.  Wind  tends  to  increase  the  tract  icni  of  tlie 
DMCA  dominated  by  completely  mixed  conditions. 

I  19.  Residence  time  distribution  data  for  DMCAs  and  similar  large, 
shallow  basins  are  scarce  and  are  of  generally  poor  quality.  The  remote  loca¬ 
tions  of  most  of  these  sites  and  the  extremely  long  sampling  times  required 
for  adequate  data  collection  are  partially  responsible.  Residence  time  dis¬ 
tribution  data  from  similar,  but  smaller  water  and  wastewater  treatment  basins 
indicate  a  strong  direct  relationship  between  basin  length-width  ratio  (h/W) 
and  hydraulic  efficiency  and  an  inverse  relationship  between  1,/W  and  the 
dispersion  index.  DMCA  hydraulic  efficiency  data  seem  to  follow  the  same 
general  relationship  with  h/W  as  data  from  the  other  sources. 

120.  The  functional  relationship  between  hvdraulic  efficiency  and  h/W 
allows  economic  optimization  of  layouts  of  DMCAs.  Use  of  one  to  three  longi¬ 
tudinal  sptir  dikes  is  economically  justified  for  most  situations.  Constrain¬ 
ing  DMCA  size  in  one  dimension  (for  example,  width)  results  Iti  higher  costs 
than  for  an  unconstrained  configuration. 

121.  Mathematical  computer  models  can  he  used  to  simulate  flow  patterns 
and  the  transport  of  solids  and  dves  in  DMCAs.  lliree  rel.ated  olistacles  t 
their  use  as  design  tools  remain,  however:  (a)  computer  costs  for  model  runs 
are  quite  high,  (b)  tlio  use  of  existing  knowledge  to  select  .appropriate  model 


coefficients  for  a  given  DMCA  is  difficult,  and  (c)  detailed,  high-quality 
field  data  are  needed  for  refinement  of  modeling  techniques. 

Recommendat ions 

122.  Part  V  of  this  report  should  be  used  when  sizing  DMCAs  until 
superseded  by  superior  information. 

123.  The  relationships  between  hydraulic  efficiency  and  L/W  developed 
in  this  document  rely  on  dye  tracer  data  from  basins  other  than  DMCAs  for 
values  of  L/W  greater  than  4.1.  Residence  time  distributions  from  DMCAs 
with  higher  L/W  should  be  measured  in  order  to  verify  the  equations  pre¬ 
sented  in  this  report.  Collection  of  highly  detailed  quality  field  data  for 
DMCAs  of  a  variety  of  shapes  would  allow  verification  of  the  relationship 
between  hydraulic  efficiency  and  L/W  ,  allow  formulation  of  design  guidance 
based  on  the  entire  residence  time  distribution  instead  of  hydraulic  effi¬ 
ciency,  and  also  allow  further  development  of  computer  modeling  capability. 

124.  If  additional  field  data  allow  formulation  of  a  relationship 
between  solids  removal  and  L/W  with  greater  precision  than  the  relationship 
between  hydraulic  efficiency  and  L/W  presented  above,  consideration  should 
be  given  to  developing  a  more  sophisticated  approach  to  DMCA  design  optimiza¬ 
tion.  An  optimization  study  might  include  consideration  of  DMCA  shape,  spur 
dikes,  design  of  spur  and  exterior  dikes,  weir  design  and  placement,  and  inlet 
structures. 

125.  Field  tests  to  determine  DMCA  residence  time  distribution  should 
include  frequent  measurements  of  wind  speed,  wind  direction,  and  flow  rate  as 
well  as  effluent  fluorescence.  Measurement  of  effluent  fluorescence  should 
commence  prior  to  dye  addition  at  the  inlet  to  allow  determination  of  back¬ 
ground  fluorescence  and  should  continue  until  fluorescence  levels  return  to 
background  levels.  Field  tests  should  also  comply  with  guidance  contained  in 
Appendix  B. 

126.  Additional  Investigations  are  needed  to  determine  the  relative 
merits  of  different  types  of  tracer  dyes  when  used  in  DMCA  tests. 

127.  Future  work  should  also  focus  on  prediction  of  DMCA  effluent 
quality  by  combining  solids  removal  curves  (Figure  1)  with  residence  time 
distributions.  Use  of  the  mean  residence  time  for  design  can  be  misleading 
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APPENDIX  A:  NUMERICAL  MODELING 


PART  I :  INTRODUCTION 


Purpose 

1.  The  purpose  of  this  study  was  to  evaluate  the  potential  of  the  TABS 
numerical  modeling  system  as  a  tool  to  predict  residence  time  distribut  ions  of 
dredged  material  containment  areas  (DMCAs) . 

Scope 

2.  The  study  consisted  of  trial  applications  of  two-dimensional  (2-D) 
and  three-dimensional  (3-n)  numerical  models  from  the  TABS  system  to  DMCAs  for 
which  measurements  of  residence  time  distribution  were  available.  This  report 
describes  the  models,  explains  how  they  were  used  to  model  DMCAs,  provides 
results  of  model  runs  for  DMCAs,  and  offers  conclusions  about  the  usefulness 
of  these  tools  in  DMCA  design. 


Approach 

3.  The  TABS  numerical  modeling  system  (Thomas  and  McAnally  1985)*  was 
used  to  simulate  eight  DMCA  test  cases.  Model  performance  was  evaluated  based 
on  two  criteria:  reproduction  of  measured  dye  tracer  curves  (which  are  resi¬ 
dence  time  d i s t r Ibu t ions)  and  prediction  of  hydraulic  efficiency.  Regression 
equations  were  derived  for  the  model  coefficients  to  allow  the  models  to  be 
used  in  a  predictive  mode.  The  models  were  used  to  simulate  a  hypothetical 
DMCA  generated  bv  adding  spur  dikes  to  one  test  case  geometry  to  show  the 
response  of  the  models  to  such  changes  in  geometrv. 


*  See  References  at  the  end  of  the  main  text. 
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PART  II:  NLT1KRICAI,  MODKI.S 


4.  Finite  element  I'.vd  rinivnamic  models,  KMA-IO  (K  ing,  198?)  and  RMA-?',' 
('Thomas  and  McAnallv  198''),  and  transport  models,  8F!)-8  (Ariathurai  19RP)  and 
RMA-4  (Thomas  and  >’cAnall\’  198S)  from  the  TABS  svstem,  were  used  in  this 
studv.  Computations  are  performed  in  two  horizontal  directions  witliin  RMA-?'.’ 
and  RMA-4.  Two-dimensional  hi'rizontal  models  are  applicable  to  I)M(,'.\s  when 
there  is  little  variation  in  velocity  in  the  vertical  plane.  Adjustment  of 
some  of  the  two-d i mens i ona 1  nuniel  coefficients  can  compensate  for  such  three- 
dimensional  phenomena  as  vertical  mixing  and  return  flow.  KM.\-1()  atid  S!?i-F 
perform  fullv  three-d  im.ens  ional  computations. 


The  Hydrodyn  amic  Models,  RMA- 1(1  and  RMA  -  ? 

5.  RMA-10  and  RM.A-i’V  are  finite  element  solutions  <tf  the  Reynolds  form 
of  the  Navier-Stokes  equ.ations  for  turbulent  tlov.  In  R.'’1A-?V,  the  equations 
are  depth-integrated.  Bottom  friction  is  calculated  with  Manning's  equation, 
and  eddy  viscosity  ct'e f f i c i ent s  are  used  to  define  the  turbulent  exchanges, 
velocity  form  of  the  basic  equation  is  used.  Side  boundaries  are  treatetl  as 
either  slip  (flow  parallel  to  the  boundary''  or  ru'-slip  (zero  flow)  as  speci¬ 
fied  by  the  user.  B(.nind.irv  conditions  mav  be  w.iter  levels,  \’e  loci  ties,  or 
discharges  and  may  occur  inside  the  mesh  as  well  tis  along  the  edges. 


he  Transp('rt  M<Hlels,  SFH-H  atid  RMA-4 


h.  The  transport  models,  SKIl-8  and  RM.A-4,  •-.olve  the  cotn.'ec  t  i  on- 
diffusion  equation,  whicli  h.as  general  :-ource-s  i  til  terms.  In  K'''A-4  ,  the  equa¬ 
tions  arc  dept  h- i  nt  eg  r.i  t  cd  ,  and  up  t<i  six  dissolved  or  susj'ended  cor,',  t  i  t  t  ' 
can  be  refuted.  In  Sld'-8  ,  t  r.insport  of  dissolved  or  suspended  raaterhil  , 
including  cohesive  sediT'ci-r,  m.r,-  be  niodeled.  ''M.'-',  us,.-  lo.  d  tod'.n.ir'  i  c  ■  ceuer- 


,ited  b\’  RM.A-?\'  .and  the  s.ur.e  coiTipr,  t  .a  t  ion.al  mesh, 
generated  bv  RMA-IH. 
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Mesh  c.eneration 

7.  For  tlie  f  i  n ;  t  (.■  element  metlinci,  the  DMCAs  were  subdivided  intn 
smaller  areas  callis!  elements.  l.'ithin  the  TAllS  svstem,  elements  can  be  either 
cnad  ri  1  ater.i  1  nr  trianjtii'ar  iir  shape  and  can  have  curved  suies.  lach  element 
is  defined  hv  a  number  c't  points  calleii  nrnies  that  are  connected  b\'  !  inee 

(element  siiies''.  'Aithin  TAHS  ,  rrnad  r  i  1  a  t  e  ra  1  s  are  defined  bv  eittht  nodes,  .in' 
trianttles  are  letined  fe.-  six  nodes.  Bed  elevations  are  def  ined  at  each  node, 
and  a  liriMr  i  n  t  i- r  pr '  I  a  t  i  im  is  performed  a  1  <'n,e,  each  element  side'.  (hitput  is 
read:!'.'  ava  i  1  a'- I  e  at  eacli  nisie,  bnt,  due  to  the  continuoirs  soliitioir  uithiii'  tin 
element,  solution  values  mat'  be  ('btained  anvwhere  within  the  ele’tent  h.s-ain. 

b,  iCr  eacli  il  i  sposa  ’  area  tested,  the  limits  of  the  area  weri-  ilicitined 
into  regions  tor  input  ti'  tlu‘  TAHS  mesh  generator  prov.ram,  AI’'!'('1‘!‘(|;  ai:d 

McAiial!''  '.s’here  allowed  'n-  area  pei'metrv,  elements  of  nnifi'rrrr  d.inaun  ioi 

!'V  s(i  !t  were  used  f  ('r  all  oi  tlie  sites  tt.odeled  in  I'fder  to  elimirrate  ele¬ 
ment  sine  as  OIK'  possil'le  source  of  a  tl  i  f  f  e  r  ence  in  eddv  visctisitv  and  di  flu- 
si  '■  coelficieiit'  ,  bed  elevations  were  specified  to  produce  a  s  1  op,-  of 
!  •' r  'l.OCK)  t' t  f  rc'm  tlie  inflow  pr'int  to  the  isitflciw  weir.  'ihe  deltaic  >■(.,]  j  [pj.;' t 

deposits  around  tlie  inf  low  pipes  were  not  modeled  ’u-canse  <d  lack  of  data 
do'SL;  r  ill  i  lu;  tlie  formati(ins.  Figure  A1  shciws  a  samt'le  ci'inpti  t  a  t  i  ona  1  mesli. 


Fieure  A!.  Numerical  model  r”<s!i,  d'est  I 
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inflow  pipe.  An  exit  water  surface  elevation  equal  to  the  DMf'A  water  surlace 
elevation  was  specified  to  simulate  the  outflow  weir.  For  the  RMA-10  test, 
slip  flow'  was  specified  at  the  subsurface  nodes  of  the  outflow  weir.  The 
hydrodvnamics  programs  were  used  to  simulate  steady-state  velocitv  and  deptlis 
fields,  which  were  then  input  to  the  transport  model. 

10.  Dye  tracer  tests  were  simulated  by  setting  an  initial  concentration 
value  of  one  at  the  inflow  and  zero  for  the  rest  of  the  DMCA.  After  one  or 
two  time-steps  (depending  on  the  testl,  the  inflow  concentration  was  reduced 
to  zero.  The  Inflow  concentration  was  left  at  one  for  enough  time  so  that 
measurable  concentrations  could  be  seen  near  the  outflow  weir.  Htuir  c'r  half- 
hour  time-steps  were  used  depending  on  the  size  of  the  DMCA. 
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PART  III;  FIPI.D  DATA 


Measured  Characteristics 

I 

[  11.  Table  A1  shows  geometric  and  hydraulic  properties  of  e.uTi  fielil 

\ 

)  test.  Test  numbers  are  the  same  as  those  used  in  the  main  bodv  o‘  t!.e  rejs  rt 

Data  from  Tests  9,  11,  12,  and  13  were  not  available  in  time  fi'r  i  luo'rpi  ■  ra  t  i  <■ 
in  this  numerical  model  study.  Methods  for  determining  site  length,  surlaee 
area,  average  depth,  discharge,  average  wind-speed,  surface  drift  vclocit'., 
and  wind  angle  are  given  in  the  main  bodv  of  this  report. 

Descriptive  Calculations 

12.  Mean  residence  time,  hydraulic  efficiencv,  modal  time,  and  dis¬ 
persion  index  were  computed  for  both  the  field  test  dve  curves  and  the  i!'’e 
curves  output  by  RMA-4  using  the  DYKCON  computer  progr.am  (Haves  Kormu 

las  .and  definitit'ns  for  these  statistical  descriptors  are  given  in  the  main 
bodv  cif  this  report  . 


PART  IV:  MODEL  RESULTS 


Three-Dimensional  Model 


13.  The  TABS-3  modeling  system  was  applied  to  Test  1.  Results  showed 
insufficient  vertical  circulation  to  warrant  the  extra  time  and  cost  of 
running  a  three-dimensional  model.  Therefore,  the  modeling  effort  was  con¬ 
tinued  using  TABS-2. 


Two-Dimensional  Model 


Model  adjustment  criteria 

14.  Coefficients  within  the  TABS  modeling  system  were  adjusted  to 
obtain  accurate  simulation  of  measured  residence  time  distributions.  The  fol¬ 
lowing  paragraphs  describe  the  adjustment  criteria  that  were  used  for  eddy 
viscosity*  in  the  hydrodynamic  model  and  diffusion  coefficients*  in  the  trans¬ 
port  model. 

15.  The  hydrodynamic  portion  of  the  modeling  system  was  adjusted  in 
order  to  produce  realistic  eddy  patterns  while  keeping  the  model  numerically 

stable.  Eddy  viscosity  coefficients  were  adjusted  in  the  hydrodynamic  model. 

2 

Values  of  0,5  Ib/sec/ft  produced  tne  most  desirable  velocity  field.  Lower 
values  caused  the  hydrodynamic  model  to  become  numerically  unstable,  whereas 
higher  values  damped  out  the  eddy  patterns.  A  sample  plot  of  the  velocity 
field  is  shown  in  Fig  ire  A2.  The  presence  of  eddies  slowed  the  simulated 
transport  of  dye, 

16.  The  transport  portion  of  the  modeling  system  was  adjusted  by  vary¬ 
ing  the  diffusion  coefficients.  Diffusion  coefficients  were  selected  for  each 
of  the  eight  field  tests  by  trial  and  error.  First,  the  diffusion  coefficient 
was  varied  to  obtain  the  best  fit  of  the  peak  of  the  simulated  curve  to  the 
observed  peak.  The  resulting  eight  values  are  hereinafter  referred  to  as  the 
best  peak  fit  (BPF)  coefficients.  Then  the  diffusion  coefficient  was  adjusted 
until  the  simulated  dye  curves  had  mean  residence  times  within  10  porcoiu  ol 
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Figure  A2 .  Velocity  vector  field,  Test  6 
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Results 


17.  The  diffusion  coefficients  that  best  fit  the  BPF  and  BHE  adjustment 
criteria  are  presented  in  Table  A2.  Table  A2  also  shows  that  when  the  model 
was  adjusted  to  meet  the  BPF  criteria,  mean  residence  times  were  predicted 
within  40  percent  of  the  observed  values  in  all  but  one  case  and  within 
14  percent  of  the  observed  value  in  half  of  the  cases.  Except  for  Test  2,  in 
cases  where  the  BPF-adjusted  model  achieved  good  fit  for  mean  residence  time, 
reasonably  good  fit  was  also  obtained  for  the  dispersion  index.  Evidently, 
the  model  provided  much  better  simulation  of  some  of  the  field  tests  than 
others,  but  the  reasons  why  are  not  obvious.  For  Tests  1,  3,  7,  and  8,  the 
BPF  and  the  BHE  adjustment  criteria  did  not  produce  the  same  diffusion  coef¬ 
ficients.  The  observed  and  simulated  dye  tracer  curves  are  shown  in 
Figure  A3.  Normalized  dye  concentrations  (the  y-coordlnate  in  Figure  A3)  were 
obtained  by  dividing  the  actual  values  by  the  peak  concentration  for  each 
curve.  Raising  the  diffusion  coefficients  moved  the  peak  of  the  dye  tracer 
curve  to  the  left  on  the  time  axis  (faster  transport),  whereas  lowering  the 
coefficients  moved  the  peak  of  the  curve  to  the  right  (slower  transport) . 


Regression  Analysis 


18.  As  noted  above,  a  common  value  for  eddy  viscosity  was  used  in  the 
hydrodynamic  model  for  all  simulations,  but  different  values  for  the  transport 
model  diffusion  coefficient  were  used  in  each  test  to  meet  the  BPF  and  BHE 
criteria  (Table  A2) .  A  multiple  linear  regression  analysis  was  performed  to 
formulate  an  equation  for  predicting  the  diffusion  coefficient  based  on 
selected  physical  characteristics  of  the  field  tests.  The  regression  analysis 
was  performed  using  the  Statistical  Package  for  the  Social  Sciences  (SPSS) 
system  of  computer  programs  (Nie  et  al.  1975). 

19.  Regression  analyses  were  performed  using  both  BPF  and  BHE  diffusion 

coefficients  as  the  dependent  variable.  A  noncoil  inear  set  of  independent 

variables,  consisting  of  the  surface  drift  velocity,  surface  area,  discharge, 

and  wind  angle,  was  selected  in  a  stepwise  manner.  Surface  area  was  selected 

instea  1  of  depth  or  length  because  it  was  uncorrelated  with  the  other  variable 

of  interest  and  because  it  was  measured  directly  from  site  pl.ins  for  all  tests 

For  the  BPF  regression,  the  independent  variables  were  added  to  the  regression 

2 

equation  one  at  a  time  until  the  coefficient  of  determination  (R  )  exceeded 
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Figure  A3.  Observed  dye  tracer  curves  (solid  lines)  and  simulated  curves  from 
TABS  model  output  (dashed  lines).  The  x-axls  is  time  after  dye  injection  in 
hours,  and  the  y-axls  is  dye  concentration  divided  by  peak  dye  concentration. 
BPF  curves  show  output  from  model  with  coefficients  adjusted  to  achieve  the 
best  fit  of  the  model  output  to  the  observed  dye  curve  peak.  BHE  curves  show 
model  output  with  coefficients  adjusted  to  obtain  values  of  hydraulic 

efficiency  close  to  observed  values 


0.950.  The  same  independent  variables  produced  the  highest  R  for  BHE  regres¬ 
sion.  Table  A3  shows  the  correlation  matrix  for  the  independent  variables 
tested . 

20.  Table  A4  shows  the  multiple  linear  regression  coefficients  resulting 
from  the  calculations.  Table  A5  shows  calculated  values  of  the  model  coeffi¬ 
cients  from  the  regression  analysis  and  the  best-fit  model  coefficients  from 


the  modeling  exercise.  Table  A6  shows  some  statistics  from  the  regression 
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analysis.  R  Is  the  coefficient  of  determination,  and  values  above  0.90 
Indicate  a  "good"  fit  to  the  data.  Lower  values  for  the  standard  deviation 
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indicate  less  error  in  fitting  the  data.  The  values  for  F  can  be  compared 
with  a  table  of  F  values  for  the  degrees  of  freedom  within  the  regression. 
For  this  regression  (4  variables,  8  tests),  F  values  above  3.84  show  a 
95-percent  confidence  level,  and  values  above  7.01  show  a  99-percent  confi¬ 
dence  level  (Miller  and  Freund  1977).  However,  reliable  regression  analysis 
requires  15  to  20  observations  per  variable.  Since  these  analyses  included  a 
total  of  5  variables,  75  to  100  observations  are  needed  to  establish  a  reli¬ 
able  regression  equation  instead  of  just  8.  Figures  A4  and  A5  show  calculated 
versus  observed  scatter  plots  for  the  best  hydraulic  efficiency  and  the  best 
curve  fit,  respectively.  The  regression  fit  was  better  using  the  BPF  diffu¬ 
sion  coefficients  and  the  resulting  regression  equation  might  be  used  as  a 
starting  point  for  selection  of  diffusion  coefficients  for  future  TABS-2 
simulations  of  DMCAs. 

BPF  ADJUSTMENT  CRITERION,  4  VARIABLE  REGRESSION 


Figure  A4 .  Calculated  versus  observed  diffusion  coefficients 
for  regression  analysis  using  BPF  adjustment 


Reasons  for  Differences  Between  Observed  and 
Simulated  Residence  Time  Distributions 


21.  Sources  of  variation  between  model  simulations  and  field  tests 


could  have  been  caused  by  improper  definition  of  the  bed  geometry  and  inflow 
or  outflow  conditions,  not  modeling  vertical  circulation,  the  presence  of 


BHE  ADJUSTMENT  CRITERION,  4  VARIABLE  REGRESSION 


Figure  A5.  Calculated  versus  observed  diffusion  coefficients 
for  regression  analysis  using  BHE  adjustment 


suspended  sediment  within  the  flow,  and  errors  in  field  data  measurement. 
Despite  the  variations,  the  TABS-2  models  were  adjusted  to  closely  reproduce 
four  of  eight  observed  dye  curves.  Specifically,  sources  of  error  included: 

a.  Inflows  for  the  larger  areas  were  unsteady  since  dredges  start 
and  stop  many  times  during  a  typical  dye  test.  The  TABS-2 
models  used  constant  inflows. 

b.  Assuming  a  uniform  slope  for  the  bed  is  an  oversimplification. 
Each  area's  bed  will  have  irregularities.  Since  sediment  is 
entering  at  the  inflow,  shoaling  of  this  area  will  occur. 

£.  Vertical  mixing  caused  by  wind  is  not  modeled  directly  within 
the  two-dimensional  models,  but  is  compensated  for  by  adjusting 
the  diffusion  coefficients. 

d.  Suspended  sediment  tends  to  make  water  denser  and  may  also 
adsorb  dye. 

e.  Some  of  the  field  dye  curves  are  noisy,  and  some  of  the 
calculations  from  them  may  not  be  meaningful.  Monitoring  of 
effluent  dye  concentration  was  not  continued  long  enough  to 
adequately  define  the  decay  limbs  of  the  curves  in  most  of  the 
tests.  Accurate  background  dye  concentration  was  not  available 
for  some  of  the  tests.  Field  data  for  Tests  7  and  8  are 
suspect  as  the  hydraulic  efficiency  was  calculated  to  be 
greater  than  100  percent. 


Geometric  Sensitivity  Test 

22.  To  further  test  the  applicability  of  the  TABS-2  modeling  system  to 
DMCA  design,  the  geometry  of  one  DMCA  was  modified  while  keeping  the  hydraulic 
properties  and  model  coefficients  constant.  Spur  dikes  were  added  to  Test  7 
(Figure  A6) .  The  dye  tracer  curve  moved  to  the  right  (slower  transport. 

Figure  A7)  when  compared  with  Test  7  without  spur  dikes  (Figure  A3) . 


Figure  A6.  Velocity  vector  field.  Test  7  with  spur  dikes 
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t  A7.  live  tracer  curve,  Test  7  with  spur  dikes 


PART  V:  CONCLUSIONS 


23.  The  tests  conducted  with  the  numerical  models  showed  the  potential 
usefulness  of  the  TABS  system  for  design  of  dredged  material  disposal  areas. 
The  following  conclusions  can  be  drawn  from  these  tests: 

a.  The  tests  clearly  showed  the  impact  of  various  eddy  viscosity 
and  diffusion  coefficients  on  the  results.  The  results  were 
sensitive  to  changes  in  both  parameters. 

b.  An  eddy  viscosity  coefficient  of  0.5  Ib/sec/ft^  produced 
reasonable  current  patterns  and  realistic  travel  times  for  the 
tracer  dye  in  the  disposal  areas  tested. 

c.  There  exists  a  definite  lower  limit  for  the  eddy  viscosity 
coefficients  where  results  produced  are  unrealistic  and 
unstable.  The  possibility  exists  that  ^or  a  given  disposal 
area  and  mesh,  a  value  of  0.5  Ib/sec/ft  may  be  too  low. 

d.  Equations  that  will  predict  the  transport  model’s  diffusion 
coefficients  have  been  derived.  The  diffusion  coefficients 
depend  upon  DMCA  surface  area,  discharge,  and  wind  conditions. 

e.  Changes  in  site  geometry  by  the  addition  of  spur  dikes  influ¬ 
ence  the  numerical  model  results  in  a  reasonable  manner. 

f.  Better  field  data  with  more  accurate  measurement  of  the  tracer 
background  concentration  and  the  decay  limb  of  the  dye  curve 
are  needed. 

g.  Ciiven  more  time,  it  is  felt  that  the  three-dimensional  TABS 
models  can  be  tuned  to  better  reproduce  the  effects  of  wind. 


APPENDIX  B;  PROCEDURES  FOR  MEASURING  DMCA  RESIDENCE  TIMES 


Fluorescent  Dyes 


General 

1.  The  most  common  method  to  measure  the  residence  time  distribution  of 
a  body  of  water  is  by  the  use  of  a  dye  tracer,  typically  a  fluorescent  dye. 
Fluorescent  materials  used  for  tracing  are  unique  in  that  they  efficiently 
convert  absorbed  light  into  emitted  light  in  a  characteristic  pair  of  spectra. 
With  the  proper  light  source  and  filters,  a  fluorometer  can  measure  small 
amounts  of  fluorescent  material  in  a  sample.  Thus,  when  a  fluorescent  dye  is 
mixed  with  a  given  parcel  of  water,  that  parcel  may  be  identified  and  traced 
through  a  water  system.  The  mean  residence  time  and  the  amount  of  mixing  of 
the  water  parcel  in  the  system  can  be  quantified  by  measuring  the  dye  concen¬ 
tration  of  the  water  leaving  the  system. 

Physical-chemical  considerations 

2.  For  a  given  fluorescent  dye,  the  interaction  of  the  dye  with  sur¬ 
rounding  environmental  conditions  should  be  considered.  Use  of  a  dye  in  fresh 
water  normally  is  not  affected  by  chemical  changes;  however,  if  the  dye  were 
to  be  used  in  waters  having  high  chloride  concentrations  (salt  water),  the  dye 
loss  could  be  significant.  Photochemical  decay  of  dye  concentration  must  also 
be  considered  when  planning  a  dye-tracer  study.  Factors  influencing  photo¬ 
chemical  decay  are  light  intensity,  cloud  cover,  water  turbidity,  and  water 
column  depth.  Other  physical-chemical  impacts  on  dyes  are  related  pH,  temper¬ 
ature,  and  salinity.  Under  acidic  conditions,  adsorption  occurs  more  freely, 
resulting  in  a  loss  of  dye  and  a  reduction  in  fluorescence.  A  general  rule  of 
thumb  on  temperature  impacts  is  that  fluorescence  of  a  given  concentration  of 
dye  decreases  5  percent  for  every  2°  C  increase  in  temperature.  Tests  have 
shown  that  dve  decay  occurs  at  a  slower  rate  under  saline  conditions  (7.02-m 
sodium  chloride  solution)  (Smart  and  Laidlaw  1977).*  Additional  gtiidance  for 
designing  dye-tracer  studies  and  details  of  physica 1 -chemi ca 1  effects  on  dyes 
are  found  in  the  following:  Abood ,  Lawler,  and  Disco  1969;  Pritchard  ami  Car¬ 
penter  I960;  Feuerstein  and  Sellcck  1963;  Watt  1965;  Smart  and  Laidlaw  1977; 
Yotsukura  and  Kilpatrick  1973;  Deaner  1973;  and  Wilson  1968. 
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3.  Fluorescent  dyes  have  been  used  since  the  early  1900s.  Several  have 
been  developed  and  used  with  varying  degrees  of  success  in  the  tracing  of  sur¬ 
face  and  ground  waters.  Sm;.  i  and  Laidlaw  (1977)  evaluated  eight  dyes: 
Fluorescein,  Rhodamine  B,  Rhodamine  WT,  Sulpho  Rhodaitiine  B,  Lissamine  FF, 
Pyramine,  Amino  G  Acid,  and  Photine  Cl'.  Rhodamine  R  is  stable  in  sunlight; 
however,  it  is  readily  adsorbed  ti'  sediments  in  water.  Rhodamine  WT  was 
developed  specifically  for  water  tracing  and  is  recommended  for  such  use  by 
Johnson  1984,  Montgomery  1978,  and  Wilson  !4hH.  Altht'ugh  Rhodamine  WT  is 
generally  thought  to  have  a  verv  low  to'icitv  level,  ingestion  should  be 
restricted  to  0.75  mg/day,  and  concent  r.it  ii'ns  should  be  limited  to  less  than 
10  ppb  near  water  Intakes  (Johnson  19H4').  Because  of  its  color,  sustained 
effluent  concentrations  above  50  ppm  should  he  avoided. 

Measurement  Techniques 


Theory  of  operation 

4.  Unlike  sophisticated  and  complex  analytical  laboratory  spectro- 
f luorometers ,  filter  fluorometers  are  relatively  simple  instruments.  Basi¬ 
cally,  filter  fluorometers  are  composed  of  six  parts:  light  (excitation 
energy)  source,  primary  or  "excitation"  filter,  sample  compartment,  secondary 
or  "emlttance"  filter,  photomultiplier,  and  readout  device. 

a.  When  a  fluorescent  material  is  placed  in  a  fluorometer,  the 
spectral  portion  of  the  light  source  that  coincides  with  the 
known  excitation  spectrum  of  the  test  material  is  allowed  to 
pass  through  the  primary  filter  to  the  sample  chamber.  This 
energy  is  absorbed  by  the  fluorescent  material ,  causing 
electrons  to  be  excited  to  higher  energy  levels.  In  returning 
to  its  ground  state,  the  fluorescent  material  emits  light  that 
is  always  at  a  longer  wavelength  and  lower  frequency  than  the 
light  that  was  absorbed.  It  is  this  property  that  is  the  basis 
of  fluorometry,  the  existence  of  a  unique  pair  of  excitation  and 
emission  spectra  for  different  fluorescent  materials.  Finally, 
only  a  certain  band  of  the  emitted  light  is  passed  through  the 
secondary  filter  to  the  photomultiplier,  where  a  readout  device 
indicates  the  relative  intensity  of  the  light  reaching  it. 

Thus,  with  different  light  sources  and  filter  combinations,  the 
fluorometer  can  discriminate  between  different  fluorescent 
materials. 

b.  The  selection  of  light  sources  and  filters  is  crucial  since  thev 
determine  the  sensitivity  and  selectivity  of  the  analysis. 
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Fluorometer  manufacturers  recommend  and  supply  lamps  and  filters 
for  most  applications,  including  Rhodamine  WT  applications. 

c.  Two  types  of  fluorometers  are  in  common  field  use  today.  The 
standard  Instrument  used  in  water  tracing  by  many  groups, 
including  the  US  Geological  Survey  (Wilson  1968)  ,  has  been  the 
Turner  Model  III  manufactured  by  G.  K.  Turner  Associates. 

Turner  Designs  has  capitalized  on  recent  advances  in  electronics 
and  optics  and  developed  a  fluorometer,  the  Model  10  series, 
that  is  better  adapted  to  field  use  than  the  Turner  Model  III. 

Field  use 

5.  Once  a  fluorometer  is  calibrated,  it  must  be  decided  where  and  how 
field  samples  will  be  analyzed — in  situ  or  in  a  laboratory,  continuously  or 
discretely.  During  in  situ  analysis,  the  operation  of  the  fluorometer  in 
flow-through  mode  (where  water  from  a  given  discharge  point  in  the  containment 
area  is  pumped  continuously  through  the  sample  chamber  in  tne  fluorometer)  is 
advantageous  over  its  operation  in  cuvette  mode  (where  a  discrete  sample  is 
analyzed).  Specifically,  in  situ  flow-through  analysis  (a)  allows  the  homo¬ 
geneity  and  variation  of  fluorescence  in  the  discharge  to  be  easily  observed 
and  (b)  eliminates  the  need  for  handling  individual  samples.  Also  during  in 
situ  flow-through  analysis,  a  strip  chart  recorder  can  be  attached  to  the 
fluorometer,  simplifying  data  collection  by  providing  a  continuous  record  of 
the  fluorescence  measured.  During  laboratory  analysis,  however,  the  flow¬ 
through  system  is  seldom  used  since  discrete  samples  are  homogeneous  and 
usually  lack  the  volume  needed  to  fill  the  system.  Instead,  the  fluorometer 
is  operated  in  cuvette  mode  where  only  a  small  portion  of  a  sample  is  required 
for  analysis.  See  Johnson  (1984)  for  a  detailed  description  of  differences 
between  models  and  operation  features. 

6.  Each  method  of  analysis  also  has  its  Inherent  problems.  laboratory 
analysis  requires  that  discrete  samples  be  collected,  bottled,  labeled,  stored 
in  the  field,  and  then  transported  to  the  laboratory;  this  introduces  many 
opportunities  for  samples  to  be  lost  through  mislabeling,  misplacement,  or 
breakage.  Also,  the  frequency  of  sampling  may  be  insufficient  to  clearly 
define  the  changes  in  dye  concentration  as  a  function  of  time. 

7.  In  situ  analysis,  on  the  other  hand,  is  usually  performed  under 
adverse  environmental  conditions — often  a  a  fast  pace,  in  a  cramped, 
makeshift  work  space,  or  in  less  than  ideal  weather  conditions.  Thus,  it  is 
more  likely  that  an  error  will  occur  during  in  situ  analysis  than  during  anal¬ 
ysis  in  the  controlled  environment  of  a  laboratory.  It  is  also  usuallv 


necessary  to  compute  and  apply  many  more  temperature  correction  factors  to 
fluorescence  values  during  in  situ  analysis  than  during  a  laboratory  analysis, 
since  the  samples  to  be  analyzed  in  situ  have  not  had  a  chance  to  reach  a  com¬ 
mon  temperature.  This  also  Increases  the  chances  for  error  during  analysis. 

In  addition,  in  situ  analysis  is  usually  final.  That  is,  if  questions  are 
raised  about  the  validity  of  a  measurement  after  the  analysis,  no  sample  is 
available  for  verification.  In  situ  analysis  may  not  be  used  when  turbidity 
interference  occurs,  as  may  be  the  case  for  the  effluent  of  many  DMCAs . 

8.  To  minimize  the  risk  involved  in  relying  on  either  method  alone,  a 
combination  of  the  two  may  be  employed — a  preliminary  in  situ  analysis  to  help 
guide  the  sampling  effort  and  a  final  laboratory  analysis  to  ensure  accurate 
results  for  quantitative  analysis. 

9.  Regardless  of  when  and  where  fluorometric  analysis  takes  place, 
general  precautionary  measures  should  he  taken  to  ensure  that  the  analysis  is 
reliable . 

a.  The  fluorometer  should  be  accurately  calibrated. 

b.  Sample  contamination  should  be  avoided  by  rinsing  or  flushing 
the  sample  chamber  between  readings. 

c.  The  fluorometer  operator  should  have  experience  with  the 
instrument  that  is  used.  Experience  can  be  gained  through 
practice  with  samples  of  known  fluorescence  prior  to  the 
analysis. 

d.  Sample  temperatures  should  be  observed  and  recorded  during 
analysis  to  determine  the  necessary  fluorescence  correction 
factors . 

e.  All  information  used  to  determine  concentration  units  should  be 
recorded  (i.e.,  scale  and  meter  or  dial  deflection). 

£.  The  calibration  should  be  checked  on  a  regular  basis  (every  hour 
or  so).  This  is  especially  important  if  the  fluorometer  is 
powered  by  a  battery;  when  the  battery  is  drained,  readings  are 
no  longer  accurate. 

10.  For  flow-through  analysis  in  particular,  all  connections  between 
the  sampling  hose,  fluorometer,  and  pump  must  be  tight  to  prevent  air  bubbles 
from  entering  the  sample  chamber.  Air  bubbles  may  also  be  introduced  by  a 
leaky  pump  seal;  thus,  it  is  recommended  that  the  pump  be  connected  to  the 
svstem  so  that  water  is  drawn  up  through  the  fluorometer  to  the  pump.  A 
screen  placed  at  the  intake  end  of  the  sampling  hose  will  prevent  sand  and 
pebbles  from  altering  the  optics  of  the  system  since  thev  mav  scratch  the 
glass  in  the  sample  chamber  as  thev  travel  tlirough  the  svstem. 


^  'V  'V  '^~^J 


11.  When  analyzing  samples  in  cuvette  mode,  the  optics  of  the  system  may 
be  distorted  by  scratches  or  smudges  on  the  cuvette,  making  it  necessary  to 
wipe  the  cuvette  clean  prior  to  its  insertion  in  the  sample  chamber.  Once  the 
cuvette  is  inside  the  warm  sample  chamber,  a  reading  must  be  made  quickly  to 
prevent  warming  of  the  sample  or  condensation  forming  on  the  cuvette.  Warming 
of  the  sample  would  cause  a  reduction  in  fluorescence,  whereas  condensation 
would  distort  the  system  optics. 


Sampling  equipment 

12.  The  basic  equipment  needed  to  perform  a  dye  tracer  study  includes 
the  following: 

a.  Fluorometers  and  accessories  (filters,  spare  lamps,  recorders, 
cuvettes,  and  sample  holders).  A  spare  fluorometer  should  be 
included  if  available  since  the  entire  field  study  centers 
around  its  operation. 

b.  Standard  dye  solutions  for  the  calibrating  fluorometers. 

c.  Generators  or  12-volt  deep-cycle  marine  batteries  (with 
charger)  to  power  fluorometers  and  pumps,  if  the  dye  concentra¬ 
tion  is  to  be  monitored  continuously. 

d.  Sampling  equipment — pump  and  hoses,  automatic  sampler  or  dis¬ 
crete  sampler  (e.g.,  a  Van  Dorn  sampler),  bottles,  labels, 
waterproof  markers. 

e.  Temperature-measuring  device  for  measuring  sample  temperatures 
if  the  temperature  of  the  samples  being  analyzed  will  vary 
significantly. 

£.  Dye,  dilution  vessels,  and  injection  equipment  (e.g.,  bucket, 
pump,  and  hoses). 

£.  Description  and  dimensions  of  the  containment  area  and  survey¬ 
ing  equipment  to  measure  dimensions  of  containment  area. 

h.  Equipment  and  records  to  determine  the  flow  rate  of  the 

effluent  from  the  containment  area  (e.g.,  production  records, 
dredge  discharge  rate,  weir  length,  depth  of  flow  over  the 
weir,  and  head  above  the  weir). 

i^.  Miscellaneous  equipment  (e.g.,  life  jackets,  tool  kits). 

2.  Data  forms. 

Additional  equipment  might  include  cameras,  radios,  rain  gear,  rope,  and 
lights.  All  equipment  should  be  checked  for  proper  performance  prior  to 
transporting  to  tiie  field. 


I 


jrrj’.jrrjv^'^ 

V 

V 


Preparatory  tasks 

13.  Prior  to  conducting  the  dye-tracer  study,  the  average  discharge  rate 
at  all  points  of  discharge  from  the  containment  area  should  be  measured  or 
estimated.  If  significant  unaccountable  leakage  is  suspected,  the  inflow  rate 
must  also  be  measured  to  determine  the  leakage  rate.  Equipment  should  be 
prepared,  calibrated,  and  Installed  to  measure  or  estimate  the  discharge  rate 
during  the  dye-tracer  study.  If  dredge  production  records  are  to  be  used  to 
estimate  the  discharge,  the  discharge  should  be  correlated  to  production.  The 
average  discharge  rate  q  is  equal  to  the  sum  of  the  average  discharge  rate 

at  each  discharge  point  including  leakage. 

14.  A  survey  of  the  containment  area  should  be  performed  to  determine 
the  area,  depth,  and  volume  of  ponding  V  at  the  site  for  determination  of 
the  theoretical  residence  time  T  .  The  volume  can  be  estimated  from  as-built 
or  design  drawings  of  the  site,  but  the  depth  of  fill  and  ponding  should  be 
verified  in  the  field  if  an  accurate  estimate  of  the  hydraulic  efficiency  is 
to  be  determined  from  the  dye-tracer  study.  The  ponded  volume  is  needed  to 
estimate  dye  requirements.  An  accurate  determination  of  the  volume  is  not 
needed  if  only  the  mean  residence  time  is  to  be  determined. 

15.  Using  the  average  discharge  rate  and  the  ponded  volume,  the  theoret¬ 
ical  residence  time  of  the  site  should  be  computed  to  plan  the  duration  of  the 
dye-tracer  study  and  to  determine  the  hydraulic  efficiency. 


T  =  V/q 


(B-1) 


16.  The  background  fluorescence  should  be  measured  at  the  site.  Back¬ 
ground  fluorescence  is  the  sum  of  all  contributions  to  fluorescence  by  mate¬ 
rials  other  than  the  fluorescent  dye.  The  best  method  to  determine  the 
background  fluorescence  is  to  measure  the  fluorescence  of  the  discharge  from 
the  site  several  times  prior  to  the  addition  of  dye  at  the  inlet.  If  the 
background  fluorescence  is  expected  to  be  variable,  the  fluorescence  of  super¬ 
natant  from  the  influent  should  be  measured  before  and  during  the  dve-tracer 
study.  The  fluorescence  of  the  water  at  the  dredging  site  should  not  be  used 
as  the  background  fluorescence  since  some  of  the  sediment  that  is  mixed  with 
the  site  water  may  remain  suspended  and  exhibit  fluorescence.  Similarly,  the 
sediment  mav  release  or  adsorb  fluorescence  materials  that  would  alter  the 
fluorescence  of  the  site  water. 


17.  The  effect  of  turbidity  on  the  measurement  of  fluorescence  should  be 
examined  to  determine  whether  the  discharge  samples  should  be  filtered  prior 

to  measuring  their  fluorescence.  Turbidity  will  reduce  the  fluorescence  by 
absorbing  and  :,catterlng  the  light  from  the  fluorometer  lamp.  Filtering  is 
necessary  only  when  samples  are  highly  turbid  or  when  the  turbidity  varies 
significantly.  The  effect  of  turbidity  can  be  tested  very  simply.  A  sample 
of  the  discharge  is  divided  in  half,  and  a  small  amount  of  dye  is  added  to  one 
of  the  portions.  The  fluorometer  is  blanked  or  zeroed  on  the  portion  without 
dye  in  it,  and  the  fluorescence  of  the  portion  containing  dye  is  measured. 

Next,  both  samples  are  filtered  or  centrifuged  to  remove  the  turbidity.  The 

process  is  then  repeated  using  the  filtrates  or  supernatants — blanking  the 
fluorometer  on  the  portion  without  dye  in  it  and  measuring  the  fluorescence  of 
the  portion  containing  dye.  If  the  measured  fluorescence  of  the  sample 

without  turbidity  differed  from  the  measured  fluorescence  of  the  sample  with 

turbidity,  then  it  is  evident  that  turbidity  affected  the  analysis.  Alterna¬ 
tively,  distilled  water  could  be  used  as  the  blank  when  the  turbidity  or  the 
background  fluorescence  is  expected  to  vary  significantly  during  the  study. 

Dye  dosage  requirements 

18.  Dye  is  usually  released  instantaneously  as  a  slug  to  measure  the 
mean  residence  time  or  hydraulic  efficiency  of  a  basin.  The  dye  marks  a  small 
parcel  of  water  that  disperses  as  the  parcel  passes  through  the  basin. 

Ideally,  the  dispersion  in  a  settling  basin  is  kept  very  low,  and  the  parcel 
moves  as  a  slug  through  the  basin  by  plug  flow.  In  practice,  the  net  flow¬ 
through  velocity  is  very  low,  sufficiently  low  that  in  the  absence  of  external 
forces  the  parcel  would  move  by  plug  flow  with  some  longitudinal  dispersion. 
However,  containment  areas  are  subiect  to  wind  forces  that  transform  the 
basins  into  partially  mixed  basins  where  the  velocities  Induced  by  wind  are 
much  greater  than  the  net  flow-through  velocity.  Consequently,  the  flow 
through  the  basin  more  closely  represents  completely  mixed  conditions  than 
plug  flow  conditions.  Therefore,  the  dye  requirements  are  determined  assuming 
that  the  dye  is  completely  mixed  in  the  basin  rather  than  longitudinally 

d  i  spersed . 

19.  A  typical  dye-tracer  curve  for  a  dredged  material  containment  area, 
shown  in  Figure  Bl,  shows  a  residence  time  distribution  that  is  characteristic 
of  a  partially  mixed  basin.  Dye  appears  qulcklv  at  the  discharge  point  at 

t Ime  t  ^  , 


and  then  shortly  thereafter  the  peak  concentration  is  discharged  at 


time  tp  .  After  the  peak  concentration  reaches  the  discharge  point,  the  dye 
concentration  quickly  decreases  to  about  30  to  60  percent  of  the  peak  concen¬ 
tration  depending  on  the  wind  and  the  theoretical  residence  time  of  the  basin. 
The  dye  concentration  then  gradually  decreases  until  all  of  the  dye  is  finally 
discharged  at  time  t^  .  The  mean  residence  time  and  theoretical  residence 
time  are  shown  in  Figure  B1  as  t  and  T  ,  respectively.  The  residence  time 
distribution  indicates  that  some  of  the  water  short-circuits  to  the  discharge 
point  before  the  dye  is  completely  mixed  throughout  the  containment  area. 
However,  the  dye  becomes  well  mixed  soon  after  the  peak  concentration  is  dis¬ 
charged,  and  then  the  dye  concentration  decreases  gradually  (instead  of 
rapidly  as  before  being  completely  mixed)  to  zero. 

20.  Before  determining  the  dye  dosage  requirements  for  a  study,  a  stan¬ 
dard  calibration  curve  should  be  developed  for  the  dye  and  fluorometers  to  be 
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Figure  Bl.  A  typical  plot  of  the  residence  time  distribution 
for  dredged  material  containment  areas 


used.  This  consists  of  plotting  the  fluorometer  response  for  at  least  five 
knovm  concentrations  of  dye.  The  design  dye  concentration  is  based  on  the 
requirement  to  produce  accurately  measurable  concentrations  of  dye  for  the 
length  of  the  study  while  not  exceeding  the  maximum  fluorometer  response  or 
excessively  coloring  the  water. 

21.  The  dye  dosage  requirements  are  based  on  achieving  an  initial  con¬ 
centration  of  30  ppb  in  a  completely  mixed  basin.  This  concentration  of 
Rhodamine  WT  corresponds  to  30  percent  of  the  full  scale  deflection  of  many 
commonly  used  f luorometers .  With  this  quantity  of  dye,  the  peak  concentration 
will  generally  be  less  than  100  ppb  (or  100  percent  of  the  maximum  fluorometer 
response)  except  for  very  small  containment  areas  (<15  acres)  or  for  areas 
with  very  bad  channeling  and  short-circuiting.  Since  the  peak  concentration 
may  exceed  the  capacity  of  the  fluorometer,  discrete  samples  should  be  taken 
during  the  period  when  the  peak  concentration  is  being  discharged.  These 
samples  may  be  diluted  to  measure  the  peak  concentration. 

22.  The  dye  dosage  requirements  are  computed  as  follows; 


Dye  Dosage,  lb  =  0.00272  (C^*  PPt>)  acre-ft)  (B-2a) 


=6.24  X  10"®  (C^,  ppb)  (Vp,  ft^)  (B-2b) 

=  2.21  X  10"^  (C^,  ppb)  Vp,  1)  (B-2c) 


where 

C  =  desired  dye  concentration  (generally  30  ppb  for  Rhodamine 
WT) 

V  =  ponded  volume 
P 

Dye  Dosage  =  quantity  in  pounds  of  pure  dye  to  be  added  to  con¬ 
tainment  area 

23.  Fluorescent  dyes  are  not  generally  produced  at  100-percent 
strength.  Rhodamine  WT  is  typically  distributed  at  20-percent  dye  by  weight. 
Consequently,  the  quantity  of  manufacturer  stock  dye  would  be  five  times  as 
large  as  computed  in  Equation  B-2. 


_ Dye  Dosage _ 

Stock  Dye  Dosage  =  Stock  Concentration 


(B-3) 


where  the  stock  concentration  is  the  fractional  dye  content  by  weight.  The 
volume  of  stock  dye  required  can  be  computed  as  follows: 


„  ,  ^  ^  Stock  Dye  Dosage 

Volume  of  Stock  Dye  =  ^  /  ,,  .  .* 

^  Specific  Weight 


Lght 


(8-4) 


The  specific  gravity  of  liquid  Rhodamlne  WT  dye  at  a  concentration  of 
20  percent  by  weight  is  about  1.19,  corresponding  to  a  specific  weight  of 
1.19  kg/£  or  9.92  Ib/gal. 

Dye  addition 

24.  The  dye  should  added  to  the  influent  stream  in  liquid  form  in  a 
quantity  and  manner  that  is  easy  to  manage.  If  the  dye  comes  in  solid  form, 
it  should  be  dissolved  prior  to  adding  it.  Solid  dye  is  easier  to  transport, 
but  it  is  often  inconvenient  to  dissolve  at  field  locations.  The  dye  may  be 
diluted  to  a  volume  that  will  ensure  good  mixing  with  the  Influent  stream,  but 
the  quantity  should  not  be  so  large  that  it  takes  more  than  about  5  or 

10  minutes  to  add  the  dye.  The  dye  may  be  pumped  into  the  influent  pipe  or 
poured  into  the  influent  jet  or  pool.  Greater  dilutions  should  be  used  to 
ensure  good  mixing  if  the  dye  is  to  be  poured  Iv.to  the  influent.  Care  must  be 
taken  that  the  dye  is  distributed  as  it  flows  into  the  containment  area  in  the 
same  manner  as  the  influent.  Caution  should  be  used  to  prevent  concentrated 
dye  from  contaminating  sampling  and  measuring  equipment. 

Sampling  procedures 

25.  Sampling  should  be  conduct>-‘d  at  all  points  of  discharge  from  the 
containment  area.  The  dye  concentration  may  be  measured  continuously  at  the 
discharge,  or  discrete  samples  may  be  collected  throughout  the  test.  Discrete 
samples  must  be  taken  when  turbidity  interference  occurs  since  the  samples 
must  be  filtered  or  centrifuged.  Discrete  samples  should  be  taken  when  the 
dye  is  being  measured  continuously  to  provide  a  backup  in  the  case  of  equip¬ 
ment  malfunction  and  to  verify  the  results  of  the  continuous  monitor. 

26.  The  sampling  frequency  should  be  scheduled  to  observe  any  sig¬ 
nificant  change  in  dye  concentration  (about  5  to  10  percent  of  the  peak  dye 
concentration) .  Sampling  should  be  more  frequent  near  the  start  of  the  test 
when  dye  starts  to  exit  from  the  containment  area  and  when  the  peak  dye  con¬ 
centration  passes  the  discharge  points.  About  40  carefully  spaced  samples 
should  clearly  define  the  residence  time  distribution  or  dye-tracer  curve. 
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27,  The  sampling  duration  should  be  sufficiently  long  to  permit  the  dye 
concentration  to  decrease  to  at  least  10  percent  of  the  peak  concentration. 
For  planning  purposes,  the  duration  should  be  two  to  three  times  the  theo¬ 
retical  residence  time. 

28.  The  flow  rate  at  all  points  of  discharge  from  th  containment  area 
should  be  measured.  If  the  flow  rate  varies  significantly  (more  than 

20  percent  of  average) ,  it  should  be  measured  periodically  throughout  the 
test.  Production  records  may  be  used  to  provide  an  indication  of  the 
variability  of  the  flow  rate.  The  rate  for  flow  over  weirs  may  be  estimated 
by  measuring  the  depth  of  flow  over  the  weir  and  the  length  of  the  weir  crest 
and  then  applying  the  weir  formula  for  sharp-crested  weirs: 


Q  =  3.3  BH- 


(B-5) 


Q  =  2.6  Bh' 


(B-6) 


where ; 


Q  =  flow  rate,  cfs 
B  =  weir  crest  length,  ft 
H  =  static  head  above  weir  crest,  ft 
h  =  depth  of  flow  above  weir  crest,  ft 


Data  analysis 


29,  Data  reduction.  The  data  should  be  tabulated  in  the  following  form 


Samole 


Time  from 
Dve  Addition 


Dye  Concentration  Time 
Above  Background  Interval 


Column  1  is  the  number  of  the  sample  1  .  If  the  dye  concentration  was  moni¬ 
tored  continuously,  discrete  points  on  the  dye  concentration  curve  may  be  used 
as  samples.  Column  2  is  the  time  t^  that  elapsed  between  the  time  that  the 
dye  was  added  to  the  influent  and  the  sample  was  taken  from  the  effluent. 
Column  3  is  the  flow  rate  at  the  time  that  the  sample  was  taken.  The 

estimate  of  flow  rate  Is  needed  only  when  this  rate  is  not  constant  during  the 


iK- 


test.  Column  4  is  the  dye  concentration  of  the  sample  discounted  for  the 
background  fluorescence  ;  that  is: 


C.  =  C  .  -  . 

1  SI  bi 


(B-7) 


where 

C.  =  dye  concentration  discounted  for  background 
fluorescence  of  sample  i 

=  measured  fluorescence  of  sample  i 

=  background  fluorescence  at  time  t^ 

If  the  background  fluorescence  does  not  vary,  C,  .  will  be  a  constant  and  may 

bi 

be  eliminated  from  the  expression  for  calculating  if  the  fluorometer  is 

blanked  or  zeroed  with  the  site  water.  Column  5  is  the  interval  of  time  At. 

1 

over  which  the  sample  is  representative  of  the  results.  The  value  of  this 
interval  is  one  half  of  the  interval  between  the  times  when  the  samples 
immediately  preceding  and  following  the  sample  of  interest  were  taken. 


At.  = 
1 


t.^1  -  t.  , 
1+1  1-1 


(B-8) 


where 

At^  =  time  interval  over  which  sample  i  is  representative 

t^^j  =  time  when  the  following  sample  was  taken  where 

t  ,  is  time  when  the  effluent  fluorescence  returns 
n+1 

to  the  background  level 

t.  ,  =  time  when  the  preceding  sample  was  taken  where  t  =0 

1-1  r  D  r  p 

A  data  table  is  produced  for  each  point  of  discharge. 

30 .  Determination  of  mean  residence  time.  After  generating  the  data 
tables,  the  mean  residence  time  is  computed  as  follows: 


^  ^ t .C.q .At . 


where : 


t  =  mean  residence  time 

n  =  total  number  of  samples  (the  sum  of  the  samples  from  each 
sampling  location) 

If  the  flow  rate  is  nearly  constant  throughout  the  test,  the  equation  may  be 
simplified  to: 


t 


n 


z 


t .C.At . 
11  1 


n 


i=0 


(B-IO) 


If  the  sampling  interval  is  constant  (i.e.,  At^  =  constant),  but  the  flow  rate 
is  not  constant,  the  equation  may  be  simplified  to: 


t 


t.C.q. 
La  1  11 

1=0 


i=0 


B-ll) 


If  both  the  sampling  interval  and  flow  rate  are  constant,  the  equation  may  be 
simplified  to: 


n 


i=0 


t  =  — 
n 


i=0 


(B-12) 


R13 


31.  An  alternative  method  of  computing  t  for  constant  flow  with  a 
variable  sampling  interval  is  presented  in  the  main  body  of  this  report.  The 
same  method  las  been  computerized  in  a  program  called  DYECON,  which  is  part  of 
a  family  of  programs  named  ADDAMS  (Automated  Dredging  and  Disposal  Alterna¬ 
tives  Management  System)  by  Hayes  et  al,  (1985). 

32.  Determination  of  hydraulic  efficiency.  The  hydraulic  efficiency  is 
the  ratio  of  the  mean  residence  time  to  the  theoretical  residence  time,  where 


Hydraulic  Efficiency  =  — 


(B-13) 


33.  The  correction  factor  for  containment  area  volume  requirements  is 
equal  to  the  inverse  of  the  hydraulic  efficiency.  This  correction  is  applied 
by  multiplying  the  volume  by  the  correction  factor. 


Hydraulic  Efficiency  Correction  = 
Factor  for  Volume  Requirements 


Hydraulic  Efficiency 


(B-14) 


34.  Other  descriptions  of  residence  time  distribution.  The  despersion 
index,  modal  time,  initial  time,  and  Morrill  index  are  other  common 
descriptions  of  residence  time  distributions.  The  DYECON  program  can  compute 
these  values.  The  equation  for  computing  the  dispersion  index  (which  may  be 
the  best  indicator  of  residence  time  distribution)  is  presented  in  the  main 
body  of  this  report. 


APPENDIX  C:  FIELD  DATA 


Field  measurements  for  the  twelve  DMCA  dye  tracer  tests  described  in  Part  III 
are  shown  in  Tables  Cl  through  C13.  The  tables  are  listings  of  data  files 
that  were  input  to  the  computer  program  DYECON  to  generate  the  statistics 
given  in  Table  Cl.  These  data  files  contain  both  measured  values  and  pre¬ 
dicted  values  for  dye  curve  decay  limbs  (see  Part  IV).  The  three  numbers  at 
the  top  of  each  listing  represent  the  number  of  data  points,  the  background 
concentration,  and  the  theoretical  residence  times  in  hours.  Dummy  theoreti¬ 
cal  residence  times  of  20  hr  were  used  for  Tests  9  and  11,  since  actual  values 
were  unknown.  Dye  concentrations  are  listed  in  the  second  column,  and  times 
are  in  the  third  column.  Dye  concentrations  are  in  units  of  fluorescence,  and 
times  are  in  hours,  except  for  Test  3,  for  which  times  are  in  minutes. 


Table  Cl 

DMCA  Dye  Tracer  Test  1 


1 

.00  51.80 

THESIS  Yi 

1 

0. 

1 .00 

2 

1.98 

6.82 

3 

3.60 

11.75 

4 

4.80 

16.45 

5 

6.00 

21.15 

6 

6.60 

23.50 

7 

7.20 

25.85 

8 

9.00 

30.55 

9 

10.20 

35.25 

10 

11.10 

40.42 

11 

12.00 

41.13 

12 

12.18 

42.54 

13 

15.00 

39.95 

14 

16.02 

41.36 

15 

16.80 

39.95 

16 

18.60 

39.48 

17 

19.20 

38.54 

18 

22.20 

28.67 

19 

24.00 

25.85 

20 

25.80 

25.38 

21 

27.00 

23.97 

22 

30.00 

24.44 

23 

31.20 

23.50 

24 

36.00 

18.10 

25 

40.20 

16.45 

26 

44.40 

16.22 

27 

69.00 

7.80 

28 

94.00 

4.20 

29 

1 18.00 

2.50 

30 

143.00 

1.50 

31 

168.00 

1  ,00 

Table  C2 


Table  C2  (Concluded) 


Table  C3 

DMCA  Dye  Tracer  Test  3 


42  0.10  3.52  BURNSVILLE  TEST  3 


1 

0. 

0.10 

2 

13.00 

0. 

3 

14.00 

1.35 

4 

16.00 

1.20 

5 

20.00 

1.24 

6 

24.00 

2.06 

7 

28.00 

2.88 

8 

32.00 

2.68 

9 

36.00 

3.72 

10 

40.00 

5.72 

1 1 

44.00 

5.46 

12 

48.00 

9.56 

13 

52.00 

21.50 

14 

54.00 

12.50 

15 

56.00 

5.80 

16 

60.00 

3.90 

I  7 

64.00 

4.06 

18 

68.00 

4.06 

19 

74.00 

4.26 

20 

80.00 

3.10 

21 

86.00 

2.20 

22 

92.00 

1 .60 

23 

98,00 

1.70 

24 

104.00 

1.50 

25 

110.00 

1  .35 

26 

118.00 

1.62 

27 

126.00 

1.80 

28 

134.00 

1.90 

29 

142.00 

1  .90 

30 

150,00 

1.80 

31 

158.00 

1.76 

32 

166.00 

1.58 

33 

174.00 

1.26 

34 

186.00 

1.20 

35 

198.00 

1  .38 

36 

210.00 

1 .16 

37 

222,00 

1 .08 

38 

264.00 

0.47 

39 

307.00 

0.31 

40 

349.00 

0.21 

41 

391 .00 

0.14 

42 

434.00 

0,10 

Table  C5 


DMCA  Dye  Trarcr  Test  5 


1  15.00 

154.94 

YAZOO  1-1 

1 

0. 

15.00 

2 

1  .00 

15.00 

3 

2.00 

35.50 

4 

3.00 

51.00 

5 

4.00 

54.50 

6 

5.00 

51.00 

7 

6.00 

67.50 

8 

7.00 

58.00 

9 

8.00 

95.00 

10 

9.00 

66.00 

1  I 

10.00 

64.50 

12 

11.00 

59.50 

13 

12.00 

62.00 

14 

13.00 

50.50 

1  5 

15.00 

59.00 

16 

17.00 

60.50 

17 

18.00 

60.00 

18 

20.00 

59.50 

19 

21.00 

60.00 

20 

23.00 

60.00 

21 

24,00 

61.50 

22 

25.00 

67.50 

23 

26.00 

60.50 

24 

28.00 

59.50 

25 

29.00 

6  5.00 

26 

30,00 

65.50 

27 

31.00 

75.00 

28 

32.00 

77.50 

29 

33.00 

89.00 

30 

34.00 

52.00 

31 

35,00 

68.00 

32 

36.00 

51.50 

33 

3  7.00 

4  5.00 

34 

38.00 

61  ,50 

35 

39.00 

68.50 

36 

40.00 

72.00 

37 

41.00 

65.50 

38 

42.00 

41.50 

39 

43.00 

70,50 

40 

45,00 

65.50 

41 

46.00 

62.50 

42 

49.00 

52.20 

43 

51.00 

52.50 

44 

52.00 

54.50 

45 

54.00 

68.50 

(Continued) 


Table  f'5  (Concluded) 


1 


77 


1.50 

225.26 

MOBILE  TEST 

1 

0. 

1.50 

2 

2.00 

1.50 

3 

4.00 

2.00 

4 

6.50 

2.00 

5 

8.00 

20.00 

6 

10.00 

1  1  .00 

7 

12.00 

46.00 

8 

16.00 

34.00 

9 

18.00 

34.00 

10 

20.00 

29.50 

1  1 

22.00 

20.50 

12 

24.00 

21  .00 

1  3 

26.00 

25.50 

14 

28.00 

23.50 

15 

30.00 

25.50 

16 

32.00 

23.50 

17 

34.00 

25.00 

18 

36.00 

24.00 

19 

38.00 

28.50 

20 

40.00 

21.00 

21 

42.00 

31  .00 

2  2 

44.00 

28.00 

2  3 

46.00 

40.50 

24 

48.00 

21  .00 

25 

50.00 

21.00 

26 

52.00 

19.00 

27 

54.00 

26.50 

28 

56.00 

19.00 

29 

58.00 

20,00 

30 

60.00 

24.50 

31 

62.00 

2  7.50 

32 

64.00 

31  ,50 

33 

66.00 

37.00 

34 

68.00 

29.00 

35 

70.00 

2  9.00 

36 

72.00 

26,00 

37 

74.00 

2  1  .  50 

38 

76.00 

3  7.50 

39 

78.00 

19.50 

40 

80.00 

19,50 

41 

8  2 . 00 

1  5 , 50 

42 

86 . 00 

1  5  ,  50 

43 

<^0.0() 

20.50 

44 

44.00 

2  7,00 

(Cont  inued’) 
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Table  C6  (Concluded) 


1 

.50  225.26 

MOBILE  TEST 

45 

98.00 

16.50 

46 

102.00 

19.50 

47 

106.00 

16.50 

48 

108.00 

15.50 

49 

112.00 

15.50 

50 

116.00 

15.00 

51 

120.00 

14.00 

52 

124.00 

10.50 

53 

132.00 

10.50 

54 

136.00 

10.50 

55 

140.00 

11.00 

56 

144.00 

14.50 

57 

148.00 

11.50 

58 

150.00 

8.50 

59 

152.00 

6.00 

60 

156.00 

10.00 

61 

160.00 

10.00 

62 

164.00 

8.50 

63 

168.00 

9.50 

64 

170.00 

8.50 

65 

172.00 

7.50 

66 

174.00 

9.50 

67 

176.00 

6.00 

68 

178.00 

10.00 

69 

184.00 

9.00 

70 

188.00 

9.50 

71 

192.00 

9.50 

72 

194.00 

8.00 

73 

254.00 

5.50 

74 

314.00 

3.00 

75 

316.00 

2.50 

76 

319.00 

2.00 

77 

320.00 

1.50 

Table  C7 

DMCA  Dye  Tracer  Test  7 


1. 00  41.41  YAZOO  1-A  //6  TEST  7 


1  0.  l.OO 

2  2.00  1.00 

3  4.00  2.00 

4  6.00  3.50 

5  8.00  5.00 

6  10.00  7.00 

7  12.00  9.20 

8  14.00  12.00 

9  16.00  16.00 

0  i:’.00  18.50 

1  18.00  17.70 

2  20.00  15.60 

3  22.00  15.20 

4  24.00  15.00 

5  26.00  14.00 

6  28.00  12.00 

7  30.00  9.50 

8  32.00  9.30 

9  34.00  9.20 

!0  36.00  9.00 

>1  38.00  8.50 

>2  40.00  8.00 

!3  42.00  7.30 

>4  44.00  6.50 

>5  67.00  5.04 

>6  89.00  3.23 

n  112.00  2.10 

28  124,20  1.72 


Table  C9 

DMCA  Dye  Tracer  Test  9 


1  .00 

20.00  N 

BLAKELY  TEST  9 

I 

0. 

1  .00 

2 

1 .50 

0. 

3 

3.75 

44.00 

4 

3.86 

73.00 

5 

4.00 

30.00 

6 

4.21 

158.00 

7 

4.38 

209.00 

8 

4.50 

200.00 

9 

4.66 

180.00 

10 

4.91 

146.00 

1 1 

5.00 

95.00 

12 

5.16 

107.00 

13 

5.33 

120.00 

14 

5.50 

120.00 

15 

5.80 

120.00 

16 

6.00 

40.00 

17 

6.16 

30.00 

18 

6.33 

16.00 

19 

6.50 

18.00 

20 

6 . 66 

12.00 

21 

7.00 

10.00 

22 

7.16 

14.00 

23 

7.33 

14.00 

24 

7.50 

14.00 

25 

7.66 

14.00 

26 

7.83 

14.00 

27 

8.00 

13.00 

28 

8.16 

14.00 

29 

8.33 

15.00 

30 

8.50 

19.00 

31 

8.66 

28.00 

32 

8.83 

25.00 

33 

9.00 

29.00 

34 

9.50 

25.00 

35 

10.00 

2.00 

36 

10.58 

6.00 

37 

11.50 

13.00 

38 

12.01 

22.00 

39 

12.58 

11.00 

40 

13.50 

17.00 

41 

14.50 

11.00 

42 

15.06 

10.00 

43 

15.55 

10.00 

44 

16.15 

9.00 

(Cent inued) 
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1.00  20.00  N  BLAKELY  TEST  9 

5  16.60  9.00 

6  17.18  8.0( 

7  17.68  8.01 

8  18.25  6.01 

9  L 

0  21 
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Table  Cll  (Concluded) 


o 

o 

20.00 

CRANEY  ISLAND  TEST  11 

45 

44.00 

13.20 

46 

45.00 

12.80 

47 

46.00 

12.00 

48 

47.00 

16.60 

49 

48.00 

11.00 

50 

49.00 

12.00 

51 

50.00 

12.00 

52 

51.00 

12.00 

53 

52.00 

13.80 

54 

53.00 

15.60 

55 

54.00 

16.70 

56 

55.00 

15.20 

57 

56.00 

13.80 

58 

57.00 

9.90 

59 

58.00 

14.50 

60 

59.00 

13.60 

61 

60.00 

12.20 

62 

62.00 

16.00 

63 

64.00 

3.70 

64 

66.00 

0.77 

65 

68.00 

0.83 

66 

70.00 

0.71 

67 

72.00 

6.70 

68 

78.00 

0.01 

69 

84.00 

2.10 

70 

90.00 

0.45 

71 

96.00 

0.21 

72 

103.00 

0.1 1 

73 

109.00 

0.07 

74 

116.00 

0.03 

Table  Cl  3 


4^1  Jt^  t»#  a-* 


DMCA  Dve  Tracer  Test 


<  4 

.00  13.52 

FWLRIV  ' 

1 

0. 

4.00 

'? 

I  .00 

4.00 

3 

2.00 

4.00 

4 

3.00 

4.00 

5 

4.00 

6.00 

6 

4.50 

57.00 

7 

5.00 

80.00 

8 

5.50 

42.00 

9 

6.00 

45.00 

10 

6.50 

67.00 

1 1 

7.00 

55.00 

12 

7.50 

37.00 

13 

8.00 

35.00 

14 

8.50 

40.00 

15 

9.00 

34.00 

16 

10.00 

19,74 

17 

I  1  .00 

12.19 

18 

12.00 

7,25 

19 

13.00 

4.20 
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